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ABSTRACT 

We report measurements of I^m, ^^a, and w from eleven supernovae at z = 0.36-0.86 with 
high-quality lightcurves measured using WFPC2 on the HST. This is an independent set of high- 
redshift supernovae that confirms previous supernova evidence for an accelerating Universe. The 
high-quality lightcurves available from photometry on WFPC2 make it possible for these eleven 
supernovae alone to provide measurements of the cosmological parameters comparable in statis- 
tical weight to the previous results. Combined with earlier Supernova Cosmology Project data, 
the new supernovae yield a measurement of the mass density ilu = 0.251:°;°^ (statistical) ±0.04 
(identified systematics), or equivalently, a cosmological constant of = 0.751qq7 (statistical) 
±0.04 (identified systematics), under the assumptions of a flat universe and that the dark en- 
ergy equation of state parameter has a constant value w = —1. When the supernova results are 
combined with independent flat-universe measurements of from CMB and galaxy redshift 
distortion data, they provide a measurement of w = — I.OS^q ^q (statistical) ±0.09 (identified 
systematic), if w is assumed to be constant in time. In addition to high-precision lightcurvc 
measurements, the new data offer greatly improved color measurements of the high-redshift su- 
pernovae, and hence improved host-galaxy extinction estimates. These extinction measurements 
show no anomalous negative E{B-V) at high redshift. The precision of the measurements is such 
that it is possible to perform a host-galaxy extinction correction directly for individual supernovae 
without any assumptions or priors on the parent E{B-V) distribution. Our cosmological fits us- 
ing full extinction corrections confirm that dark energy is required with P{^a > 0) > 0.99, a 
result consistent with previous and current supernova analyses which rely upon the identification 
of a low-extinction subset or prior assumptions concerning the intrinsic extinction distribution. 
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1. Introduction 
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Five years ago, the Supernova Cosmology 
Project (SCP) and the High-Z Supernova Search 
Team both presented studies of distant Type la 
supernovae (SNe la) in a series of reports, which 
gave strong evidence for an acceleration of the Uni- 
verse's expansion, and hence for a non-zero cosmo- 
logical constant, or dark energy density (Pcrlmut- 
ter et al. 1998; Garnavich et al. 1998a; Schmidt 
et al. 1998; Riess et al. 1998; Perlmutter et al. 
1999; for a review, see Perlmutter & Schmidt 
2003). These results ruled out a flat, matter- 
dominated (r^M = 1, r^A = 0) universe. For a 
flat universe, motivated by inflation theory, these 
studies yielded a value for the cosmological con- 
stant of Oa — 0.7. Even in the absence of as- 
sumptions about the geometry of the Universe, the 
supernova measurements indicate the existence of 
dark energy with greater than 99% confidence. 

The supernova results combined with obser- 
vations of the power spectrum of the Cosmic 
Microwave Background (CMB) (e.g., Jaffe et al. 
2001; Bennett et al. 2003; Spergel et al. 2003), 
the properties of massive clusters (e.g.. Turner 
2001; Allen, Schmidt, & Fabian 2002; BahcaU et 
al. 2003), and dynamical redshift-space distortions 
(Hawkins et al. 2002) yield a consistent picture of a 
flat universe with f^M — 0.3 and f^A — 0.7 (Bahcall 
et al. 1999). Each of these measurements is sen- 
sitive to different combinations of the parameters, 
and hence they complement each other. More- 
over, because there are three diff'erent measure- 
ments of two parameters, the combination pro- 
vides an important consistency check. While the 
current observations of galaxy clusters and dy- 
namics, and of high-redshift supernovae, primarily 
probe the "recent" Universe at redshifts of z < 1, 
the CMB measurements probe the early Universe 
at z ^ 1100. That consistent results are obtained 
by measurements of vastly different epochs of the 
Universe's history suggests a vindication of the 
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standard model of the expanding Universe. 

In the redshift range around z = 0.4-0.7, the 
supernova results are most sensitive to a linear 
combination of and f2A close to J7m — ^^A- In 
contrast, galaxy clustering and dynamics are sensi- 
tive primarily to alone, while the CMB is most 
sensitive to Om + r^A- Although combinations of 
other measurements lead to a separate confirma- 
tion of the Universe's acceleration (e.g., Efstathiou 
et al. 2002), taken alone it is the supernovae that 
provide the best direct evidence for dark energy. 
Therefore, it is of importance to improve the pre- 
cision of the supernova result, to confirm the re- 
sult with additional independent high-redshift su- 
pernovae, and also to limit the possible effects of 
systematic errors. 

Perlmutter et al. (1997, 1999) and Ricss et 
al. (1998) presented extensive accounts of, and 
bounds for, possible systematic uncertainties in 
the supernova measurements. One obvious pos- 
sible source of systematic uncertainty is the effect 
of host-galaxy dust. For a given mass density, the 
effect of a cosmological constant on the magni- 
tudes of high-redshift supernovae is to make their 
observed brightnesses dimmer than would have 
been the case with $1a = 0. Dust extinction from 
within the host galaxy of the high-redshift super- 
novae could have a similar effect; however, normal 
dust will also redden the colors of the supernovae. 
Therefore, a measurement of the color of the high- 
redshift supernovae, compared to the known col- 
ors of low-redshift SNc la. has been usc^d to pro- 
vide an upper limit on the effect of host-galaxy 
dust extinction, or a direct measurement of that 
extinction which may then be corrected. Uncer- 
tainties on extinction corrections based on these 
color measurements usually dominate the statisti- 
cal error of photometric measurements. Previous 
analyses have either selected a low-extinction sub- 
set of both low- and high-redshift supernovae and 
not applied corrections directly ("Fit C," the pri- 
mary analysis of P99), or have used an asymmetric 
Bayesian prior on the intrinsic extinction distribu- 
tion to limit the propagated uncertainties from er- 
rors in color measurements (Riess et al. 1998; "Fit 
E" of P99). 

In Sullivan et al. (2003), we set stronger limits 
on the effects of host-galaxy extinction by com- 
paring the extinction, cosmological parameters, 
and supernova peak magnitude dispersion for sub- 



sets of the SCP supernovae observed in different 
types of host galaxies, as identified from both HST 
imaging and Keck spectroscopy of the hosts. We 
found that supernovae in early-type (E and SO) 
galaxies show a smaller dispersion in peak magni- 
tude at high redshift, as had previously been seen 
at low redshift (e.g. Wang, Hoeflich, & Wheeler 
1997). This subset of the P99 sample — in hosts 
unlikely to be strongly affected by extinction — 
independently provided evidence at the 5(t level 
that > in a flat Universe and conflrmed that 
host-galaxy dust extinction was unlikely to be a 
significant systematic in the results of P99, as had 
been suggested previously (e.g.. Rowan- Robinson 
2002). The natural next step following the work 
of Sullivan et al. (2003) presented in the current 
paper — is to provide high-quality individual unbi- 
ased E{B-V) measurements that allow us to di- 
rectly measure the effect of host-galaxy extinction 
on each supernova event without resorting to a 
prior on the color excess distribution. 

The current paper presents eleven new super- 
novae discovered and observed by the SCP at red- 
shifts 0.36 < z < 0.86, a range very similar to 
that of the 42 high-redshift supernovae reported 
in Perlmutter et al. (1999; hereafter P99). The su- 
pernovae of that paper, with one exception, were 
observed entirely with ground-based telescopes; 11 
of the 14 new supernovae reported by Riess et al. 
(1998) were also observed from the ground. The 
eleven supernovae of this work have lightcurves in 
both the R and I bands measured with the Wide- 
Field/Planetary Camera (WFPC2) on the Hubble 
Space Telescope (HST), and represent the largest 
sample to date of HST-measured SNe la at high 
redshift. 

The HST provides two primary advantages for 
photometry of point sources such as supernovae. 
First the sky background is much lower, allowing 
a much higher signal-to-noise ratio in a single ex- 
posure. Second, because the telescope is not lim- 
ited by atmospheric seeing, it has very high spa- 
tial resolution. This helps the signal-to-noise ratio 
by greatly reducing the area of background emis- 
sion which contributes to the noise of the source 
measurement, and moreover simplifies the task of 
separating the variable supernova signal from the 
host galaxy. With these advantages, the preci- 
sion of the lightcurve and color measurements is 
much greater for the eleven supernovae in this pa- 
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per than was possible for previous ground-based 
observations. These eleven supernovae themselves 
provide a high-precision new set of high-redshift 
supernovae to test the accelerating universe re- 
sults. Moreover, the higher precision lightcurve 
measurements in both R- and /-bands allow us 
to make high-quality, unbiased, individual host- 
galaxy extinction corrections to each supernova 
event. 

We first describe the PSF-fit photometry 
method used for extracting the lightcurves from 
the WFPC2 images (§ 2.1). Next, in § 2.2, we 
describe the lightcurve fitting procedure, includ- 
ing the methods used for calculating accurate K- 
corrections. So that all supernovae may be treated 
consistently, in § 2.3 we apply the slightly updated 
iV'-correction procedure to all of the supernovae 
used in P99. In § 2.4, the cosmological fit method- 
ology we use is described. In § 3, we discuss the 
evidence for host-galaxy extinction (only signifi- 
cant for three of the eleven new supernovae) from 
the R-I lightcurve colors. In § 4.1, we present the 
measurements of the cosmological parameters 
and f^A from the new dataset alone as well as com- 
bining this set with the data of P99. In § 4.2, we 
perform a combined fit with our data and the high- 
redshift SNe of Riess et al. (1998). Finally, in § 4.3 
we present measurements of w, the dark energy 
equation of state parameter, from these data, and 
from these data combined with recent CMB and 
galaxy redshift distortion measurements. These 
discussions of our primary results are followed by 
updated analyses of systematic uncertainties for 
these measurements in § 5. 

2. Observations, Data Reduction, and 
Analysis 

2.1. WFPC2 Photometry 

The supernovae discussed in this paper are 
listed in Table 1. They were discovered during 
three different supernova searches, following the 
techniques described in Perlmutter et al. (1995, 
1997, 1999). Two of the searches were con- 
ducted with the 4m Blanco telescope at the Cerro 
Tololo Inter- American Observatory (CTIO), in 
November/December 1997 and March/April 1998. 
The final search was conducted at the Canada- 
France-Hawaii Telescope (CFHT) on Mauna Kea 
in Hawaii in April/May 2000. In each case, 2-3 



nights of reference images were followed 3-4 weeks 
later by 2-3 nights of search images. The two im- 
ages of each search field were seeing-matched and 
subtracted, and were searched for residuals indi- 
cating a supernova candidate. Weather conditions 
limited the depth and hence the redshift range of 
the March/ April 1998 search. Out of the three 
searches, eleven of the resulting supernova discov- 
eries were followed with extensive HST photome- 
try. These supernovae are spaced approximately 
evenly in the redshift range 0.3 < z < 0.9. Nine 
out of the eleven supernovae were discovered very 
close to maximum light; two were discovered sev- 
eral days before maximum light. 

Spectra were obtained with the red side of LRIS 
on the Keck 10m telescope (Oke et al. 1995), with 
FORSl on Antu (VLT-UTl) (Appenzeller et al. 
1998), and with EFOSC224 on the ESO 3.6m tele- 
scope. These spectra were used to confirm the 
identification of the candidates as SNe la, and 
to measure the redshift of each candidate. Nine 
of the eleven supernovae in the set have strong 
confirmation as Type la through the presence of 
Si II A6150, Si II A4190, or Fe II features that 
match those of a Type la observed at a similar 
epoch. SNe 1998ay and 1998be have spectra which 
are consistent with SNe la spectra, although this 
identification is less secure for those two. How- 
ever, we note that the colors (measured at multiple 
epochs with the HST lightcurves) are inconsistent 
with other non-la types. (We explore the system- 
atic effect of removing those two supernovae from 
the set in § 5.2.) 

Where possible, the redshift, z, of each candi- 
date was mcasiired by matching narrow features 
in the host galaxy of the supernovae; the precision 
of these measurements in z is typically 0.001. In 
cases where there were not sufficient host-galaxy 
features (SNcl998aw and 1998ba), redshifts were 
measured from the supernova itself; in these cases, 
z is measured with a (conservative) precision of 
0.01 (Branch & van den Bergh 1993). Even in 
the latter case, redshift measurements do not con- 
tribute significantly to the uncertainties in the 
final cosmological measurements since these are 
dominated by the photometric uncertainties. 

Each of these supernovae was imaged with 
two broadband filters using the Planetary Cam- 

'http://www.ls.eso.org/lasilla/sciops/efosc/ 
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Table 1: WFPC2 Supernova Observations 



"SW i 

Name Observations Observations 

1997ek OMS 1998-01-05 (400s,400s) 1998-01-05 (500s,700s) 

1998-01-11 (400s,400s) 1998-01-11 (500s,700s) 



1998-02-02 (1100s, 1200s) 
1998-02-14 (1100s, 1200s) 
1998-02-27 (1100s, 1200s) 
1998-11-09 (1100s,1300s) 
1998-11-16 (1100s,1300s) 



1997eq 0.538 1998-01-06 (300s. 300s) 1998-01-06 (300s, 300s) 

' 's.400s) 1998-01-11 (300s, 300s) 

1998-02-02 (500s. 700s) 
's,400s) 1998-02-11 (500s, 700s) 

is,400s) 1998-02-19 (500s,700s) 



1998- 


-01- 


06 


1998- 


-01- 


■21 


1998- 


-02- 


11 


1998- 


■02- 


■19 



1997ez 0.778 1998-01-05 (400s,400s) 1998-01-05 (500s, 700s) 

1998-01-11 (400s,400s) 1998-01-11 (500s, 700s) 



1998-02-02 (1100s. 1200s) 
1998-02-14 (1100s. 1200s) 
1998-02-27 (100s, 1200s, 1100s, 1200s) 



1998-04-08 (400s,400s) 

1998-04-20 (400s. 400s) 
1998-05-11 (400s. 400s) 
1998-05-15 (400s, 400s) 
1998-05-29 (400s,400s) 



1998-04-08 (500s,700s) 

1998-04-20 (500s, 700s) 

1998-05-11 (500s, 700s) 

1998-05-15 (500s, 700s) 

1998-05-29 (500s,700s) 



1998-04-08 (300s, 300s) 
1998-04-18 (300s. 300s) 
1998-04-29 (400s. 400s) 
1998-05-14 (400s. 400s) 
1998-05-28 (400s. 400s) 



1998-04-08 (300s, 300s) 
1998-04-18 (300s, 300s) 
1998-04-29 (500s, 700s) 
1998-05-14 (500s, 700s) 
1998-05-28 (500s, 700s) 



1998-04-08 (300s. 300s) 
1998-04-18 (300s. 300s) 
1998-04-29 (300s, 300s) 
1998-05-14 (300s, 300s) 
1998-05-27 (300s,300s) 



1998-04-08 (300s, 300s) 
1998-04-18 (300s, 300s) 
1998-04-29 (500s, 700s) 
1998-05-14 (500s,700s) 
1998-05-27 (500s,700s) 



1998ay 0.638 1998-04-08 (400s,400s) 1998-04-08 (500s,700s) 

1998-04-20 (400s,400s) 1998-04-20 (500s,700s) 

1998-05-11 (1100s, 1200s) 
1998-05-15 (1100s. 1200s) 
1998-06-03 (1100s. 1200s) 



1998-04-08 (300s. 300s) 

1998-04-19 (300s, 300s) 

1998-04-29 (400s, 400s) 

1998-05-13 (400s, 400s) 

1998-05-28 (400s,400s) 



1998-04-08 (300s, 300s) 
1998-04-19 (300s, 300s) 
1998-04-29 (500s, 700s) 
1998-05-13 (500s,700s) 
1998-05-28 (500s,700s) 



1998-04-08 (300s,300s) 
1998-04-19 (300s, 300s) 
1998-04-30 (400s. 400s) 
1998-05-15 (400s. 400s) 
1998-05-28 (400s. 400s) 



1998-04-08 (300s, 300s) 

1998-04-19 (300s, 300s) 

1998-04-30 (500s, 700s) 

1998-05-15 (500s, 700s) 

1998-05-28 (500s, 700s) 



1998bi 0.740 1998-04-06 (400s,400s) 1998-04-06 (500s, 700s) 

1998-04-18 (400s,400s) 1998-04-18 (500s,700s) 

1998-04-28 (1100s,1200s) 
1998-05-12 (1100s,1200s) 
1998-06-02 (1100s,1200s) 



2000-05-15 (600s, 600s) 
2000-05-28 (600s, 600s) 
2000-06-10 (500s, 500s) 
2000-06-22 (1100s,1300s) 
2000-07-08 (1100s,1300s) 



2000-05-08 (2200s) 
2000-05-15 (1100s, 1100s) 
2000-05-28 (600s, 600s) 
2000-06-10 (600s, 600s) 
2000-06-22 (1100s, 1200s) 
2000-07-08 (110s,1200s) 
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era (PC) CCD of the WFPC2 on the HST, which 
has a scale of 0.046"/pixel. Table 1 lists the dates 
of these observations. The F675W and F814W 
broadband filters were chosen to have maximum 
sensitivity to these faint objects, while being as 
close a match as practical to the rest-frame B and 
V filters at the targeted rcdshifts. (Note that all 
of our WFPC2 observing parameters except the 
exact target coordinates were fixed prior to the su- 
pernova discoveries.) The effective system trans- 
mission curves provided by STScI indicate that, 
when used with WFPC2, F675W is most similar 
to ground-based R band while F814W is most sim- 
ilar to ground-based / band. These filters roughly 
correspond to redshifted B- and V-band filters for 
the supernovac at z < 0.7, and redshifted U- and 
B- band filters for the supernovae at z > 0.7. 

The HST images were reduced through the 
standard HST "On-The-Fly Reprocessing" data 
reduction pipeline provided by the Space Tele- 
scope Science Institute. Images were then back- 
ground subtracted, and images taken in the same 
orbit were combined to reject cosmic rays using the 
"crrej" procedure (a part of the STSDAS IRAF 
package). Photometric fluxes were extracted from 
the final images using a PSF-fitting procedure. 
Traditional PSF fitting procedures assume a sin- 
gle isolated point source above a constant back- 
ground. In this case, the point source was super- 
imposed on the image of the host galaxy. In all 
cases, the supernova image was separated from the 
core of the host galaxy; however, in most cases the 
separation was not enough that an annular mea- 
surement of the background would be accurate. 
Because the host-galaxy flux is the same in all of 
the images, we used a PSF fitting procedure that 
fits a PSF simultaneously to every image of a given 
supernova observed through a given photometric 
filter. The model we fit was: 

fi{x, y) = foi X psf (a; - xoi,y - yoi) + 

hg{x-xoi,y-yoi;aj)+pi (1) 

where fi{x,y) is the measured flux in pixel (x,y) 
of the ith image, (xoi,yoi) is the position of the 
supernova on the ith image, /oi is the total flux 
in the supernova in the ith image, psf(u,w) is a 
normalized point spread function, bg(w, v\ a) is a 
temporally constant background parametrized by 
aj, and pi is a pedestal offset for the ith image. 



There arc 4n + m — 1 parameters in this model, 
where n is the number of images (typically 2, 5, or 
6 previously summed images) and m is the num- 
ber of parameters aj that specify the background 
model (typically 3 or 6). (The —1 is due to the 
fact that a zeroth-order term in the background is 
degenerate with one of the pi terms.) Parameters 
varied include /j, XQi, yoi, Pi, and aj. 

Due to the scarcity of objects in our PC images, 
geometric transformations between the images at 
different epochs using other objects on the four 
chips of WFPC2 together allowed an a priori de- 
termination of {xoi.yoi} good to ~ 1 pixel. Allow- 
ing those parameters to vary in the fit (effectively, 
using the point source signature of the supernova 
to determine the offset of the image) provided po- 
sition measurements a factor of ^ 10 better. The 
model was fit to 13 x 13 pixel patches extracted 
from all of the images of a time sequence of a single 
supernova in a single filter (except for SN 1998ay, 
which is close enough to the host galaxy that a 
7x7 pixel patch was used to avoid having to fit 
the core of the galaxy with the background model) . 
In four out of the 99 patches used in the fits to 
the 22 lightcurves, a single bad pixel was masked 
from the fit. The series of /o, values, corrected 
as described in the rest of this section, provided 
the data used in the lightcurve fits described in 
§ 2.2. For one supernova (SN 1997ek at z = 0.86), 
the F814W background was further constrained 
by a supernova-free "final reference" image taken 
11 months after the supernova explosion. 

A single Tiny Tim PSF was used as psf (m, v) 
for all images of a given filter. The Tiny Tim PSF 
used was subsampled to 10 x 10 subpixels; in the 
fit procedure, it was shifted and integrated (prop- 
erly summing fractional subpixels). After shifting 
and resampling to the PC pixel scale, it was con- 



'Note that this may introduce a bias towards higher fiux, 
as the fit wiU seek out positive fluctuations on which to 
center the PSF. However, the covariancc between the peak 
flux and position is typicaUy less than ~ 4% of the product 
of the positional uncertainty and the flux uncertainty, so 
the effects of this bias will be very small in comparison to 
our photometric errors. 

Although obtaining final references to subtract the galaxy 
background is standard procedure for ground-based pho- 
tometry of high-redshift supernovae, the higher resolution 
of WFPC2 provides sufficient separation between the su- 
pernova and host galaxy that such images are not always 
necessary, particularly in this redshift range. 
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volvcd with an empirical 3x3 electron diffusion 
kernel with 75% of the flux in the central element 
(Pruchter 2000).^'^ The PSF was normalized in a 
0.5"-radius aperture, chosen to match the stan- 
dard zeropoint calibration (Holtzman, et al. 1995; 
Dolphin 2000). Although the use of a single PSF 
for every image is an approximation — the PSF of 
WFPC2 depends on the epoch of the observation 
as well as the position on the CCD — this approxi- 
mation should be valid, especially given that for all 
of the observations the supernova was positioned 
close to the center of the PC. To verify that this 
approximation is valid, we reran the PSF fitting 
procedure with individually generated PSFs for 
most supernovae; we also explored using a super- 
nova spectrum instead of a standard star spectrum 
in generating the PSF. The measured fluxes were 
not significantly diff'erent, showing differences in 
both directions generally within 1 2% of the su- 
pernova peak flux value — much less than our pho- 
tometric uncertainties on individual data points. 

Although one of the great advantages of the 
Hubble Space Telescope is its low background, 
CCD photometry of faint objects over a low back- 
ground suffer from an imperfect charge transfer 
efficiency (CTE) effect, which can lead to a sys- 
tematic underestimate of the flux of point sources 
(Whitemore, Heyer, & Casertano 1999: Dolphin 
2000, 2003). On the PC, these effects can be 
as large as ^ 15%. The measured flux val- 
ues {foi above) were corrected for the CTE of 
WFPC2 following the standard procedure of Dol- 
phin (2000).^^ Uncertainties on the CTE correc- 
tions were propagated into the corrected super- 
nova fluxes, although in all cases these uncertain- 
ties were smaller than the uncertainties in the raw 
measured flux values. Because the host galaxy is a 
smooth background underneath the point source, 
it was considered as a contribution to the back- 
ground in the CTE correction. For an image which 
was a combination of several separate exposures 
within the same orbit or orbits, the CTE calcula- 
tion was performed assuming that each SN image 
had a measured SN flux whose fraction of the total 
flux was equal to the fraction of that individual im- 

Scc also http:/ /www. stsci.edu/softwarc/tinytim/tinytim_faq. html 

'These CTE corrections used updated co- 
efficients posted on Dolphin's web page 
(htt p : // www .noao.edu/staiT/dolphin/ wfpc2_calib/) in 
September, 2002. 



age's exposure time to the summed image's total 
exposure time. This assumption is correct most of 
the time, with the exception of the few instances 
where Earthshine affects part of an orbit. 

In addition to the HST data, there exists 
ground-based photometry for each of these super- 
novae. This includes the images from the search it- 
self, as well as a limited amount of follow-up. The 
details of which supernovae were observed with 
which telescopes are given with the lightcurves in 
Appendix A. Ground-based photometric fltixes 
were extracted from images using the same aper- 
ture photometry procedure of P99. A complete 
lightcurve in a given filter {R or /) combined 
the HST data with the ground-based data (us- 
ing the color-correction procedure described be- 
low in § 2.3), using measured zeropoints for the 
ground-based data and the Vega zeropoints of Dol- 
phin (2000) for the HST data. The uncertainties 
on those zeropoints (0.003 for F814W or 0.006 
for F675W) were added as correlated errors be- 
tween all HST data points when combining with 
the ground-based lightcurve. Similarly, the mea- 
sured uncertainty in the ground-based zeropoint 
was added as a correlated error to all ground- 
based fluxes. Ground-based photometric calibra- 
tions were based on observations of Landolt (1992) 
standard stars observed on the same photomcit- 
ric night as a supernova observation; each cal- 
ibration is confirmed over two or more nights. 
Ground-based zeropoint uncertainties are gener- 
ally < 0.02-0.03; the ii-band ground based zero- 
point for SN 1998ay is only good to ±0.05. We 
have compared our ground-based aperture pho- 
tometry with our HST PSF-fitting photometry us- 
ing the limited number of sufficiently bright stars 
present in the PC across the eleven SNe fields. We 
find the difference between the HST and ground- 
based photometry to be 0.02 ± 0.02 in both the R- 
and I-bands, consistent with no offset. The cor- 
related uncertainties between different supernovae 
arising from ground-based zeropoints based on the 
same calibration data, and between the HST su- 
pernovae (which all share the same zeropoint), 
were included in the covariance matrix used in all 
cosmological fits (see § 2.4). 

2.2. Lightcurve Fits 

It is the magnitude of the supernova at its 
lightcurve peak that serves as a "calibrated can- 
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die" in estimating the cosmological parameters 
from the luminosity distance relationship. To esti- 
mate this peak magnitude, we performed template 
fits to the time scries of photometric data for each 
supernova. In addition to the eleven supernovae 
described here, lightcurve fits were also performed 
to the supernovae from P99, including 18 super- 
novae from Hamuy et al. (1996a; hereafter H96), 
and eight from Riess et al. (1999a; hereafter R99) 
which match the same selection criteria used for 
the H96 supernovae (having data within six days 
of maximum light and located at cz > 4000 km/s, 
limiting distance modulus error due to peculiar 
velocities to less than 0.15 magnitudes). Because 
of new templates and ii'-corrections (see below), 
lightcurve fits to the supernovae from H96 and P99 
used in the analyses below were redone for consis- 
tency. The results of these fits are slightly different 
from those quoted in P99 for the same supernovae 
as a result of the change in the lightcurve template, 
the new ii'-corrections, and the different fit proce- 
dure, all discussed below. For example, because 
the measured E{B-V) value was considered in the 
iiT-corrections (§ 2.3), whereas it was not in P99, 
one should expect to sec randomly distributed dif- 
ferences in fit supernova lightcurve parameters due 
to scatter in the color measurements. 

Lightcurve fits were performed using a x^' 
minimization procedure based on MINUIT (James 
& Roos 1975). For both high- and low-redshift 
supernovae, color corrections and -fT-corrections 
are applied (see § 2.3) to the photometric data. 
These data were then fit to lightcurve templates. 
Fits were performed to the combined R- and I- 
band data for each high-redshift supernova. For 
low-redshift supernovae, fits were performed using 
only the B- and T^-band data (which correspond 
to de-redshifted R- and /-bands for most of the 
high-redshift supernovae). The lightcurve model 
fit to the supernova has four parameters to mod- 
ify the lightcurve templates: time of rest-frame B- 
band maximum light, peak flux in R, R-I color at 
the epoch of rest-frame maximum _B-band light, 
and timescale stretch s. Stretch is a parameter 
which linearly scales the time axis, so that a su- 
pernova with a high stretch has a relatively slow 
decay from maximum, and a supernova with a low 
stretch has a relatively fast decay from maximum 
(Pcrlmutter et al. 1997; Goldhaber et al. 2001). 
For supernovae in the redshift range z = 0.3-0.7, 



a B template was fit to the _R-band lightcurve and 
a V template was fit to the /-band lightcurve. For 
supernovae at ^; > 0.7, a U template was fit to 
the /?-band lightcurve and a B template to the 
/-band lightcurve. Two of the high-redshift su- 
pernovae from P99 fall at 2; ~ 0.18 (SN 19971 and 
SN1997N); for these supernovae, V and R tem- 
plates were fit to the R- and /-band data. (The 
peak B-band magnitude was extracted by adding 
the intrinsic SN la B-V color to the fit V"-band 
magnitude at the epoch of B maximum.) 

The B template used in the lightcurve fits was 
that of Goldhaber et al. (2001). For this paper, 
new F-band and i?-band templates were generated 
following a procedure similar to that of Goldhaber 
et al. (2001), by fitting a smooth parametrized 
curve through the low-redshift supernova data of 
H96 and R99. A new JZ-band template was gen- 
erated with data from Hamuy et al. (1991), Lira 
et al. (1998), Richmond et al. (1995), Suntzeff et 
al. (1999), and Wells et al. (1994); comparison of 
our 17-band template shows good agreement with 
the new JJ-band photometry from Jha (2002) at 
the relevant epochs. New templates were gener- 
ated by fitting a smooth curve, f{t'), to the low- 
redshift lightcurve data, where t' = t/ {l+z)/ s; t is 
the number of observer-frame days relative to the 
epoch of the i?-baiid maximum of each supernova, 
z is the redshift of each supernova, and s is the 
stretch of each supernova as measured from the 
B-band lightcurves. Lightcurve templates had an 
initial parabola with a 20-day rise time (Aldering, 
Knop, & Nugent 2000), joined to a smooth spline 
section to describe the main part of the lightcurve, 
then joined to an exponential decay to describe 
the final tail at >~ 70 days past maximum light. 
The first 100 days of each of the three templates 
is listed in Table 2. 

Due to a secondary "hump" or "shoulder" ^ 20 
days after maximum, the i?-band lightcurve does 
not vary strictly according to the simple time-axis 
scaling parametrized by stretch which is so suc- 
cessful in describing the different U-, B-, and V- 
band lightcurves. However, for the two z ^ 0.18 
supernova to which we fit an i?-band template, 
the peak R- and /- band magnitudes are well con- 
strained, and the stretch is also well measured 
from the rest-frame F-band lightcurve. 

Some of the high-redshift supernovae from P99 
lack a supernova-free host-galaxy image. These 
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Table 2: IJ . \ \ and R Liglitciirvo Toinplatcs Used 



Day'' 


U flux'' 


V flux'' 


R flux'= 


7^ 1 

Dayi 




U flux'' 


V flux'' 




R flux'' 


-19 


6.712e-03 


4.960e-03 


5.779e-03 


31 


4, 


.790e-02 


2.627e-01 


3 


.437e-01 


-18 


2.685e-02 


1.984e-02 


2.312e-02 


32 


4, 


.5246-02 


2.481e-01 


3 


.238e-01 


-17 


6.041e-02 


4.464e-02 


5.201e-02 


33 


4, 


.300e-02 


2.345e-01 


3 


.054e-01 


-16 


1.074e-01 


7.935e-02 


9.246e-02 


34 


4 


.1120-02 


2.2I80-OI 


2 


.887e-01 


-15 


1.678O-01 


I.24O0-OI 


1. 4450-01 


35 


3 


.9560-02 


2.099O-01 


2 


.7330-01 


-14 


2. 4160-01 


1.785O-01 


2.O8O0-OI 


36 


3 


.8270-02 


1.990C-01 


2 


.592e-01 


-13 


3.289O-01 


2.43O0-OI 


2.832O-01 


37 


3 


.7220-02 


1.891O-01 


2 


.4630-01 


-12 


4.296e-01 


3.174e-01 


3.698e-01 


38 


3, 


.636e-02 


1.802e-01 


2 


.345e-01 


-11 


5.437e-01 


4.017e-01 


4.681e-01 


39 


3 


.565e-02 


1.721e-01 


2, 


.237e-01 


-10 


6.712e-01 


4.960e-01 


5.779e-01 


40 


3 


.506e-02 


1.649e-01 


2 


.137e-01 


-9 


7.486e-01 


5.889e-01 


6.500e-01 


41 


3 


.456e-02 


1.583e-01 


2 


.046e-01 


-8 


8.151e-01 


6.726e-01 


7.148e-01 


42 


3 


.410e-02 


1.524e-01 


1 


.962e-01 


-7 


8.711C-01 


7.469e-01 


7.725e-01 


43 


3, 


.365e-02 


1.471e-01 


1 


.884e-01 


-6 


9.168C-01 


8. 1150-01 


8.236e-01 


44 


3 


.3I80-O2 


1.423O-01 


1 


.8130-01 


-5 


9.524C-01 


8.66O0-OI 


8.681e-01 


45 


3, 


.2660-02 


1.378C-01 


1 


.7470-01 


-4 


9.781C-01 


9. 1030-01 


9.062e-01 


46 


3, 


.205e-02 


1.337C-01 


1 


.6870-01 


-3 


9.940e-01 


9.449e-01 


9.382e-01 


47 


3, 


.139e-02 


1.299e-01 


1, 


.630e-01 


-2 


l.OOOe-FOO 


9.706e-01 


9.639e-01 


48 


3 


.072e-02 


1.263e-01 


1 


.578e-01 


-1 


9.960e-01 


9.880e-01 


9.834e-01 


49 


3 


.005e-02 


1.229e-01 


1 


.529e-01 





9.817e-01 


9.976e-01 


9.957e-01 


50 


2 


.945e-02 


1.195e-01 


1 


.483e-01 


1 


9.569e-01 


l.OOOe-l-00 


l.OOOe-l-00 


51 


2 


.893e-02 


1.161e-01 


1 


.440e-01 


2 


9.213e-01 


9.958e-01 


9.952e-01 


52 


2, 


.853e-02 


1.128e-01 


1 


.398e-01 


3 


8.742C-01 


9.856O-01 


9. 8030-01 


53 


2 


.8300-02 


1.096C-01 


1 


.359e-01 


4 


8.172C-01 


9. 7020-01 


9.545C-01 


54 


2 


.8270-02 


1.064O-01 


1 


.320e-01 


5 


7.575C-01 


9. 5020-01 


9.196C-01 


55 


2 


.8490-02 


1.033O-01 


1 


.2820-01 


6 


6.974e-01 


9.263e-01 


8.778e-01 


56 


2 


.793e-02 


1.003e-01 


1 


.244e-01 


7 


6.375e-01 


8.991e-01 


8.313e-01 


57 


2 


.738e-02 


9.743e-02 


1 


.207e-01 


8 


5.783e-01 


8.691e-01 


7.821e-01 


58 


2 


.684e-02 


9.467e-02 


1 


.170e-01 


9 


5.205e-01 


8.369e-01 


7.324e-01 


59 


2, 


.630e-02 


9.207e-02 


1 


.133e-01 


10 


4.646e-01 


8.031e-01 


6.842e-01 


60 


2, 


.578e-02 


8.964e-02 


1, 


.097e-01 


11 


4.113e-01 


7.683e-01 


6. 3960-01 


61 


2, 


.527e-02 


8.74I0-O2 


1 


.O6I0-OI 


12 


3.610e-01 


7.33O0-OI 


6.007O-01 


62 


2, 


.4770-02 


8. 5380-02 


1 


.0260-01 


13 


3.145O-01 


6. 9770-01 


5.691O-01 


63 


2 


.4280-02 


8.3590-02 


9 


.9IO0-O2 


14 


2.725C-01 


6.629O-01 


5. 4440-01 


64 


2 


.3800-02 


8.207O-02 


9 


.5680-02 


15 


2.356e-01 


6.293e-01 


5.254e-01 


65 


2, 


.333e-02 


8.083e-02 


9, 


.232e-02 


16 


2.044e-01 


5.972e-01 


5.113e-01 


66 


2, 


.287e-02 


7.927e-02 


8, 


.902e-02 


17 


1.783e-01 


5.667e-01 


5.011e-01 


67 


2, 


.242e-02 


7.774e-02 


8, 


.579e-02 


18 


1.567O-01 


5.376O-01 


4.938e-01 


68 


2 


.1970-02 


7.624O-02 


8 


.2640-02 


19 


1.388C-01 


5. 0990-01 


4.887e-01 


69 


2 


.1540-02 


7. 4760-02 


7. 


.9580-02 


20 


1.239C-01 


4.835O-01 


4.848O-01 


70 


2, 


.III0-O2 


7. 3320-02 


7. 


.6600-02 


21 


1.115e-01 


4.583O-01 


4.814O-01 


71 


2, 


.0700-02 


7.I9I0-O2 


7. 


.3730-02 


22 


I.OO80-OI 


4. 3420-01 


4.776O-01 


72 


2 


.0290-02 


7. 0520-02 


7. 


.0960-02 


23 


9.144C-02 


4. 1130-01 


4.725O-01 


73 


1 


.9890-02 


6.9I60-O2 


6 


.8320-02 


24 


8.314e-02 


3.894e-01 


4.653e-01 


74 


1 


.949e-02 


6.782e-02 


6, 


.581e-02 


25 


7.583e-02 


3.685e-01 


4.552e-01 


75 


1 


.911e-02 


6.651e-02 


6, 


.344e-02 


26 


6.941e-02 


3.486e-01 


4.414e-01 


76 


1, 


.873e-02 


6.523e-02 


6 


.199e-02 


27 


6. 3800-02 


3.296O-01 


4.247O-01 


77 


1 


.836e-02 


6. 3970-02 


6 


.0570-02 


28 


5.891O-02 


3. 1150-01 


4.058O-01 


78 


1 


.7990-02 


6.274O-02 


5. 


.9180-02 


29 


5.467C-02 


2. 9430-01 


3.855O-01 


79 


1 


.7640-02 


6.153O-02 


5 


.7830-02 


30 


5.102C-02 


2.781O-01 


3.645O-01 


80 


1 


.7290-02 


6.034O-02 


5 


.6500-02 



a: Day is relative to the epoch of the maximum of the S-band lightcurve. The S-band template may 
be found in Goldhaber et al. (2001). 
b: Relative fluxes. 
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supernovac were fit with an additional variable pa- 
rameter: the zero-level of the I-band lightcurve. 
The supernovae treated in this manner include 
SNe 19970, 1997Q, 1997R, and 1997am. 

The late-time lightcurve behavior may bias the 
result of a lightcurve fit (Aldering, Knop, & Nu- 
gent 2000); it is therefore important that the low- 
and high-redshift supernovae are treated in as con- 
sistent a manner as possible. Few or none of 
the high-redshift supernovae have high-precision 
measurements more than ^40-50 rest-frame days 
after maximum light, so as in Perlmutter et al. 
(1997) and P99 these latc-timc points were elim- 
inated from the low-redshift lightcurve data be- 
fore the template-fit procedure. Additionally, to 
allow for systematic offset uncertainties on the 
host-galaxy subtraction, an "error floor" of 0.007 
times the maximum lightcurve flux was applied; 
any lightcurve point with an uncertainty below 
the error floor had its uncertainty replaced by that 
value (Goldhaber et al. 2001). 

The final results of the lightcurve fits, including 
the effect of color corrections and ii'-corrections, 
are listed in Table 3 for the eleven supernovae 
of this paper. Table 4 shows the results of new 
lightcurve fits to the high-redshift supernovae of 
P99 used in this paper (see § 2.5), and Table 5 
shows the results of lightcurve fits for the low- 
redshift supernovae from H96 and R99.^^ Ap- 
pendix A tabulates all of the lightcurve data for 
the eleven HST supernovae in this paper. The 
lightcurves for these supernovae (and the F675W 
WFPC2 image nearest maximum light) are shown 
in Figures 1 and 2. Note that there are correlated 
errors between all of the ground-based points for 
each supernova in these figures, as a single ground- 
based zeropoint was used to scale each of them 
together with the HST photometry. 

2.3. Color- and ii'- Correct ions 

In order to combine data from different tele- 
scopes, color corrections were applied to remove 
the differences in the spectral responses of the fil- 
ters relative to the Bessell system (Bessell 1990). 
For the ground-based telescopes, the filters are 
close enough to the standard Bessell filters that a 
single linear color term (measured at each observa- 

' These three tables are available in electronic form from 
http: / / supernova.lbl.gov. 



tory with standard stars) suffices to put the data 
onto the Bessell system, with most corrections be- 
ing smaller than 0.01 magnitudes. The WFPC2 
filters arc different enough from the ground-based 
filters, however, that a linear term is not suffi- 
cient. Moreover, the differences between a SN la 
and standard star spectral energy distribution are 
significant. In this case, color corrections were cal- 
culated by integrating template SN la spectra (de- 
scribed below) through the system response. 

In order to perform lightcurve template fit- 
ting, a cross-filter K-correction must be applied 
to transform the data in the observed filter into 
a rest-frame magnitude in the filter used for the 
lightcurve template (Kim, Goobar, & Perlmutter 
1996). The color correction to the nearest stan- 
dard Bessell filter followed by a iiT-correction to 
a rest-frame filter is equivalent to a direct K- 
correction from the observed filter to the stan- 
dard rest-frame filter. In practice, we perform 
the two steps separately so that all photometry 
may be combined to provide a lightcurve effec- 
tively observed through a standard (e.g. i?-band) 
filter, which may then be fit with a single se- 
ries of if-corrections. The data tabulated in Ap- 
pendix A have all been color-corrected to the stan- 
dard Bessell filters. 

Color and X-corrcctions were performed follow- 
ing the procedure of Nugent, Kim, & Perlmutter 
(2002). In order to perform these corrections, a 
template SN la spectrum for each epoch of the 
lightcurve, as described in that paper, is neces- 
sary. The spectral template used in this present 
work began with the template of that paper. To 
it was applied a smooth multiplicative function at 
each day such that integration of the spectrum 
through the standard filters would produce the 
proper intrinsic colors for a Type la supernova (in- 
cluding a mild dependence of those intrinsic colors 
on stretch). 

The proper intrinsic colors for the supernova 
spectral template were determined in the BVRI 
spectral range by smooth fits to the low-redshift 
supernova data of H96 and R99. For each color 
{B-V, V-R, and R-I), every data point from those 
papers was i^-corrected and corrected for Galactic 
extinction. These data were plotted together, and 
then a smooth curve was fit to the plot of color 
versus date relative to maximum. This curve is 
given by two parameters, each of which is a func- 
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Fig. 1. Lightcurvcs and images from the PC CCD on WFPC2 for the HST supernovac reported in this 
paper. The left column shows the i?-band (including F675W HST data), and the middle column shows I- 
band lightcurves (including F814W HST data). Open circles represent ground-based data points, and filled 
circles represent WFPC2 data points. Note that there arc correlated errors between all of the ground-based 
points for each supernova in these figures, as a single ground-based zeropoint was used to scale each of them 
together with the HST photometry. The right column shows 6" x 6" images, summed from all HST images 
of the supernova in the indicated filter. 
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Fig. 2. — Lightcurves and images from tlic PC CCD on WFPC2 for the HST supcrnovac reported in this 
paper (continued). The left column shows the i?-band (including F675W HST data), and the middle column 
shows /-band lightcurves (including F814W HST data). Open circles represent groimd-based data points, 
and filled circles represent WFPC2 data points. Note that there are correlated errors between all of the 
ground-based points for each supernova in these figures, as a single ground-based zeropoint was used to scale 
each of them together with the HST photometry. The right column shows 6" x 6" images, summed from all 
HST images of the supernova in the indicated filter. 
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Table 3: Supernova Lightcurve Fits: HST Supernovae from this paper 



SN 




z 


mx 
(a) 




m_B 
(b) 






(c) 




Ext 


. Corr. 


(d) 


Stretch (s) 




R-I 

(e) 




E(B-V) 
Gal. (/) 


(9) 


>st 


Excluded from 
Subsets {h) 


1997ek 


0, 


.863 


23.32 


24. 


.51 ± 0, 


.03 


24, 


.59 ± 0. 


19 


24. 


,95 ± 0, 


.44 


1.056 ± 0.058 


0, 


.838 ± 0, 


.054 


0.042 


-0.091 ± 0, 


.075 




1997eq 


0, 


.538 


22.63 


23. 


.21 ± 0, 


.02 


23, 


.15 ± 0. 


,18 


23. 


,02 ± 0, 


.17 


0.960 ± 0.027 


0, 


.202 ± 0, 


.030 


0.044 


0.035 ± 0, 


.034 




1997ez 


0, 


.778 


23.17 


24. 


.29 ± 0, 


.03 


24, 


.41 ± 0. 


,18 


24. 


,00 ± 0, 


.42 


1.078 ± 0.030 


0, 


.701 ± 0, 


.048 


0.026 


0.095 ± 0, 


.068 




1998as 


0, 


.355 


22.18 


22. 


.72 ± 0, 


.03 


22, 


.66 ± 0. 


.17 


22. 


.02 ±0, 


.15 


0.956 ±0.012 


0, 


.226 ± 0, 


.027 


0.037 


0.158 ±0. 


.030 


2,3 


1998aw 





.440 


22.56 


23. 


.22 ± 0, 


.02 


23 


.26 ± 0. 


.17 








1.026 ± 0.019 


0, 


.300 ± 


.024 


0.026 


0.259 ±0, 


.026 


1-3 


1998ax 


0, 


.497 


22.63 


23, 


,25 ± 0, 


.05 


23, 


.47 ± 0, 


,17 


22, 


,96 ± 0, 


.20 


1.150 ± 0.032 


0, 


.212 ± 0, 


.041 


0.035 


0.113 ± 0, 


.044 


2,3 


1998ay 


0, 


.638 


23.26 


23, 


,86 ± 0, 


.08 


23, 


.92 ± 0, 


,19 


23, 


,85 ± 0, 


.33 


1.040 ± 0.041 


0, 


.339 ± 0, 


.067 


0.035 


0.015 ± 0, 


.084 


3 


1998ba 


0, 


.430 


22.34 


22, 


97 ± 0, 


.05 


22 


.90 ± 0, 


,18 


22, 


,75 ± 0, 


.18 


0.954 ± 0.020 


0, 


,094 ± 0, 


.036 


0.024 


0.040 ± 0, 


.038 




1998bo 


0, 


.644 


23.33 


23, 


91 ± 0, 


.04 


23, 


.64 ± 0, 


,18 


23, 


,26 ± 0, 


.27 


0.816 ± 0.028 


0, 


,436 ± 0, 


.051 


0.029 


0.106 ± 0, 


.065 


3 


1998bi 


0, 


.740 


22.86 


23, 


,92 ± 0, 


.02 


23, 


.85 ± 0, 


,17 


23, 


,75 ± 0, 


.37 


0.950 ± 0.027 


0, 


,552 ± 0, 


.037 


0.026 


0.026 ± 0, 


.050 




2000fr 


0, 


.543 


22.44 


23, 


07 ± 0, 


.02 


23, 


.16 ± 0, 


,17 


23, 


,27 ± 0, 


.14 


1.064 ± 0.011 


0, 


,135 ± 0, 


.022 


0.030 


-0.031 ± 


.025 





a: Magnitude in the observed filter at the peak of the rest-frame B-band lightcurve. X— i? for z < 0.7, X—/ for z > 0.7. 

b: This value has been iC-corrected and corrected for Galactic extinction: ms = 'mx — ^BX — Ax, where Kbx is the cross-filter X-correction and Ax is the Galactic 
extinction correction. These were the values used in the cosmological fits. The quoted error bar is the uncertainty on the peak magnitude from the lightcurve fit. 
c: This value includes the stretch correction: m^J^ = rriB + Ci{s — 1). a is the best-fit value of the stretch-luminosity slope from the fit to the primary low-extinction subset 
(Fit 3 in § 4). The quoted error bar includes all uncertainties for non-extinction-corrected fits described in § 2.4. Note that these values are only provided for convenience; 
they were not used directly in any cosmological fits, since a is also a fit parameter. 

d: Similar to column c, only with the host-galaxy extinction correction applied. The stretch/luminosity slope used for this value is that from the fit to the primary subset 
(Fit 6 in § 4). The quoted error bar includes all uncertainties for extinction-corrected fits described in § 2.4. A line indicates a supernova which did not appear in the primary 

subset (sec § 2.5.) 

e: This is the observed R-I color at the epoch of the rest-frame S-band lightcurve peak. 

/: Schlegel, Finkbeiner, & Davis (1998); this extinction is already included in the quoted values of ms- 

g: Measurement uncertainty only; no intrinsic color dispersion included. 

h: These supernovae are excluded from the indicated subsets; see § 2.5. 



Table 4: Supernova Lightcurve Fits: New Fits to Perlmutter (1999) SNe 



SN 


z 


mx 

(a) 




ms 

(b) 






(c) 




Ext 


m'J' 
. Corr. 


id) 


Stretch (s) 




R-I 

(e) 




E{B-V) 
Gal. (/) 


(a) 


Excluded from 
Subsets ill) 


1995ar 


0.465 


22.80 


23, 


.48 ± 0, 


,08 


23, 


,35 ± 0, 


22 


21 


,54 ± 0, 


,97 


0.909 ± 0, 


104 


0, 


,509 ± 0, 


,222 





,022 


0, 


,448 ± 0.242 




1995as 


0.498 


23.03 


23, 


.69 ± 0, 


,07 


23, 


,74 ± 0, 


,23 


23, 


,52 ± 


.87 


1.035 ± 0, 


,090 


0, 


,155 ± 0, 


,197 





,021 


0, 


.051 ± 0.212 


3 


1995aw 


0.400 


21.78 


22, 


.28 ± 0, 


.03 


22, 


.57 ±0. 


,18 


23, 


.17 ±0, 


.45 


1.194 ± 0. 


.037 


-0, 


.127 ±0, 


.103 


0, 


.040 


— 0. 


.160 ± 0.107 




1995ax 


0.615 


22.56 


23 


.21 ± 0, 


.06 


23, 


.38 ± 0. 


,22 


23 


.98 ± 1, 


.02 


1.112 ± 0. 


.073 


0, 


.152 ± 0, 


.204 


0, 


.033 


— 0. 


.153 ± 0.249 




1995ay 


0.480 


22.64 


23, 


.07 ±0, 


.04 


22, 


,90 ± 0. 


19 


22, 


.74 ± 0, 


.70 


0.880 ± 0. 


,064 


0, 


,209 ± 0, 


.158 


0, 


.114 


0. 


.047 ± 0.170 




1995az 


0.450 


22.46 


22, 


.70 ± 0, 


,07 


22, 


,66 ± 0. 


,20 


23, 


.04 ± 


.58 


0.973 ± 0. 


,064 


0, 


,087 ±0, 


,135 


0, 


.181 


— 0. 


,089 ± 0.144 




1995ba 


0.388 


22.07 


22, 


.64 ± 0, 


,06 


22, 


,60 ± 0, 


18 


22, 


,74 ± 


.45 


0.971 ± 0, 


,047 


0, 


,006 ± 0, 


,105 





,018 


~0, 


,033 ± 0.110 




1996cf 


0.570 


22.71 


23, 


.31 ± 0, 


,03 


23, 


,30 ± 0, 


18 


23, 


,53 ± 


.45 


0.996 ± 0, 


,045 


0, 


,162 ± 0, 


,091 





,040 


— 0, 


,054 ± 0.107 


3 


1996cg 


0.490 


22.46 


23, 


.09 ± 0, 


,03 


23, 


,11 ± 0, 


18 


22 


,26 ± 


.45 


1 nil -t- n 


040 


0, 


,300 ± 0, 


,099 


n 
u 


,uoo 


Q 


9nr^ -t- n 1 n? 


3 


1996ci 


0.495 


22.19 


22 


.83 ± 0, 


,02 


22, 


,78 ± 0, 


18 


22 


,92 ± 


.32 


0.964 ± 0, 


,040 


0, 


,083 ± 0, 


,070 


0, 


,028 


-0, 


,033 ± 0.075 




1996cl 


0.828 


23.37 


24, 


.53 ± 0, 


,17 


24, 


,49 ± 0, 


,46 


25, 


,92 ± 


.97 


0.974 ± 0, 


,239 


0, 


,549 ± 0, 


,184 


0, 


,035 


-0, 


,344 ± 0.251 




1996cm 


0.450 


22.67 


23, 


.26 ± 0, 


.07 


23, 


.11 ± 0. 


,18 


22, 


.63 ± 0, 


.77 


0.899 ± 0, 


.061 


0, 


.214 ±0, 


.174 


0, 


.049 


0, 


.124 ±0.185 


3 


1996cn 


0.430 


22.58 


23 


.25 ± 0, 


,03 


23, 


,09 ± 0, 


19 








0.890 ± 0, 


,066 


0, 


,379 ± 0, 


,090 


0, 


.025 


0, 


,332 ± 0.097 


1-3 


1997F 


0.580 


22.93 


23, 


.51 ± 0, 


,06 


23, 


,57 ± 0, 


,20 


23, 


,30 ± 


.95 


1.041 ± 0, 


,066 


0, 


,275 ± 0, 


,197 


0, 


,040 


0, 


,063 ± 0.232 




1997H 


0.526 


22.70 


23, 


.26 ± 0, 


,04 


23, 


,09 ± 0, 


19 


22 


,51 ± 


.80 


0.882 ± 0, 


,043 


0, 


,303 ± 0, 


,174 


0, 


,051 


0, 


,150 ± 0.194 




19971 


0.172 


20.18 


20, 


.34 ± 0, 


,01 


20, 


,29 ± 0, 


17 


20, 


.19 ± 0, 


.28 


0.967 ± 0, 


,009 


0, 


065 ± 0, 


,047 


0, 


,051 


0, 


026 ± 0.064 




1997N 


0.180 


20.39 


20, 


.38 ± 0, 


,02 


20, 


,48 ± 0, 


17 


21, 


.28 ± 0, 


.52 


1.067 ± 0, 


,015 


-0, 


,141 ± 0, 


093 


0, 


,031 


-0, 


,200 ± 0.123 




19970 


0.374 


22.99 


23, 


.53 ± 0, 


,06 


23, 


,60 ± 0, 


18 








1.048 ± 0, 


,054 


0, 


087 ± 0, 


152 


0, 


,029 


0, 


049 ± 0.162 


1-3 


1997P 


0.472 


22.53 


23, 


.16 ± 0, 


,04 


22, 


,99 ± 0, 


18 


23, 


,24 ± 


,91 


0.888 ± 0, 


039 


0, 


058 ± 0, 


207 


0, 


,033 


-0, 


052 ± 0.219 




1997Q 


0.430 


22.01 


22, 


.61 ± 0, 


,02 


22, 


,52 ± 0, 


17 


22, 


,55 ± 


.62 


0.935 ± 0, 


,024 


0, 


,061 ± 0, 


,140 


0, 


,030 


-0, 


,002 ± 0.148 




1997R 


0.657 


23.29 


23, 


.89 ± 0, 


.05 


23, 


.80 ±0. 


19 


23, 


.68 ± 0, 


.90 


0.940 ± 0, 


.059 


0, 


.393 ± 0, 


.175 


0, 


.030 


0, 


.032 ± 0.222 




1997ac 


0.320 


21.42 


21, 


.87 ±0, 


,02 


21, 


,96 ± 0. 


,17 


21, 


.95 ± 0, 


.33 


1.061 ± 0, 


,015 


0, 


,063 ± 0, 


,065 


0, 


.027 


0, 


,001 ± 0.072 




1997af 


0.579 


22.94 


23, 


.60 ± 0, 


,07 


23, 


,38 ± 0, 


18 


24, 


,31 ± 1 


.09 


0.850 ± 0, 


,045 


0, 


,045 ± 0, 


,226 


0, 


,028 


-0, 


,215 ± 0.265 




1997ai 


0.450 


22.34 


22 


.94 ± 0, 


,05 


22, 


,63 ± 0, 


,22 


22 


,58 ± 


.59 


0.788 ± 0, 


,084 


0, 


,143 ± 0, 


,133 


0, 


,045 


0, 


,026 ± 0.142 




1997aj 


0.581 


22.58 


23, 


.24 ± 0, 


,07 


23, 


,16 ± 0, 


18 


24, 


,05 ± 


.79 


0.947 ± 0, 


,045 


0, 


045 ± 0, 


,164 


0, 


,033 


-0, 


,213 ± 0.193 




1997am 


0.416 


22.01 


22 


,58 ± 0, 


,08 


22, 


,63 ± 0, 


18 


22 


,65 ± 


.46 


1.032 ± 0, 


,060 


0, 


037 ± 0, 


,113 


0, 


,036 


-0, 


,008 ± 0.119 




1997ap 


0.830 


23.16 


24, 


.35 ± 0, 


,07 


24, 


,38 ± 0, 


18 


23, 


,74 ± 


,50 


1.023 ± 0, 


,045 


0, 


,903 ± 0, 


,082 


0, 


,026 


0, 


,155 ± 0.118 





a: X^R for z < 0.7, X^I for z > 0.7 

b: This value has been i^-corrected and corrected for Galactic extinction: tub = rnx — i^BX — Ax, where Kbx is the cross-filter X-correction and Ax is the Galactic 
extinction correction. These were the values used in the cosmological fits. The quoted error bar is the uncertainty on the peak magnitude from the lightcurve fit. 
c: This value includes the stretch correction: m^J^ = ms ot{s — 1). a is the best-fit value of the stretch-luminosity slope from the fit to the primary low-extinction subset 
(Fit 3 in § 4). The quoted error bar includes all uncertainties for non-extinction-corrected fits described in § 2.4. Note that these values are only provided for convenience; 
they were not used directly in any cosmological fits, since a is also a fit parameter. 

d: Similar to column c, only with the host-galaxy extinction correction applied. The stretch/luminosity slope used for this value is that from the fit to the primary subset 
(Fit 6 in § 4). The quoted error bar includes all uncertainties for extinction-corrected fits described in § 2.4. A line indicates a supernova which did not appear in the primary 

subset (sec § 2.5.) 

e: This is the observed R-I color at the epoch of the rest-frame S-band lightcurve peak. 

/: Schlegel, Finkbeiner, & Davis (1998); this extinction is already included in the quoted values of ms- 

g: Measurement uncertainty only; no intrinsic color dispersion included. 

h: These supernovae are excluded from the indicated subsets; see § 2.5. 



Table 5: Supernova Lightcurve Fits: Low-z SNe from Hamuy (1996) and Riess (1999) 



SN 


z 




niB 














Stretch (s) 




R,-I 


E(B-V) 






Excluded from 


(a) 




(b) 


(c) 






(d) 




Ext, corr. 


(e) 






if) 


Gal. (g) 




{"■) 




Subsets (i) 


1990O 


0.030 


16.58 


16.18 ± 0, 


.03 


16, 


,33 ± 0, 


20 


16.30 ± 0, 


17 


1.106 ± 0.026 


0, 


,043 ± 0,025 


0.098 





,001 ± 


,026 




1990af 


0.050 


17.92 


17.76 ± 0, 


.01 


17, 


,39 ± 0, 


,18 


17.42 ± 0, 


,13 


0.749 ± 0.010 


0, 


,077 ± 0.011 


0.035 





,011 ± 


,011 




1992P 


0.026 


16.12 


16.05 ± 0, 


.02 


16, 


,14 ± 0, 


,19 


16.16 ± 0, 


,16 


1.061 ± 0.027 


-0, 


,045 ±0.018 


0.020 


— 0. 


,008 ± 0, 


,019 




1992ae 


0.075 


18.59 


18.42 ± 0, 


.04 


18, 


,35 ± 0, 


,18 


18.35 ± 0, 


,15 


0.957 ±0.018 


0, 


,098 ± 0.028 


0.036 


0, 


,003 ± 0, 


,031 




1992ag 


0.026 


16.67 


16.26 ± 0, 


.02 


16, 


,34 ± 0, 


,20 


15.55 ± 0, 


,16 


1.053 ± 0.015 


0, 


,220 ± 0.020 


0.097 


0, 


,189 ± 0, 


,021 


2,3 


1992al 


0.014 


14.61 


14.48 ± 0, 


.01 


14, 


,42 ± 0, 


,23 


14.53 ± 0, 


,20 


0.959 ± 0.011 


-0, 


,054 ± 0.012 


0.034 


— 


,025 ± 


,013 




1992aq 


0.101 


19.38 


19.30 ± 0, 


.02 


19, 


,12 ± 0, 


,17 


19.24 ± 0, 


,15 


0.878 ± 0.017 


0, 


,142 ± 0.023 


0.012 


— 


,019 ± 


,026 




1992bc 


0.020 


15.18 


15.10 ± 0, 


.01 


15, 


,18 ± 0, 


,20 


15.36 ± 0, 


,16 


1,053 ± 0,006 


-0, 


,087 ± 0.009 


0.022 


— 


,046 ± 


,009 




1992bg 


0.036 


17.41 


16.66 ± 0, 


.04 


16, 


,66 ± 0, 


,20 


16.68 ± 0, 


,16 


1,003 ± 0,014 


0, 


,128 ± 0.025 


0.181 


— 


,006 ± 


,026 




1992bh 


0.045 


17.71 


17.60 ± 0, 


,02 


17, 


,64 ± 0, 


,18 


17.22 ± 0, 


,14 


1,027 ± 0,016 


0, 


,101 ± 0.018 


0.022 


0, 


,100 ± 0, 


,019 




1992bl 


0.043 


17.37 


17.31 ± 0, 


.03 


17, 


,03 ± 0, 


,18 


17.10 ± 0, 


,14 


0.812 ± 0.012 


0, 


,017 ± 0.023 


0.012 


-0, 


,002 ± 0, 


,024 




1992bo 


0.018 


15.89 


15.78 ± 0, 


.01 


15, 


,42 ± 0, 


,21 


15.31 ± 0, 


,17 


0.756 ± 0.005 


0, 


,048 ± 0.012 


0.027 


0, 


,043 ± 0, 


,012 




1992bp 


0.079 


18.59 


18.29 ± 0, 


.01 


18, 


,16 ± 0, 


,18 


18.41 ± 0, 


,13 


0.906 ± 0.014 


0, 


,088 ± 0.015 


0.068 


-0, 


,056 ± 0, 


,017 




1992br 


0.088 


19.52 


19.37 ± 0, 


.08 


18, 


,93 ± 0, 


,20 






0.700 ± 0.021 


0, 


,186 ± 0.047 


0.027 


0, 


,030 ± 0, 


,052 


1-3 


1992bs 


0.063 


18.26 


18.20 ± 0, 


.04 


18, 


,26 ± 0, 


,18 


18.37 ±0, 


,14 


1.038 ± 0.016 


0, 


,011 ± 0.022 


0.013 


-0, 


,031 ± 0, 


,024 




1993B 


0.071 


18.74 


18.37 ± 0, 


,04 


18, 


,40 ± 0, 


,18 


18.10 ± 0, 


,15 


1.021 ± 0.019 


0, 


,181 ± 0.027 


0.080 


0, 


,071 ± 0, 


,029 




19930 


0.052 


17.87 


17.64 ± 0, 


,01 


17, 


53 ± 0, 


,18 


17.61 ± 0, 


,13 


0.926 ± 0.007 


0, 


,042 ± 0.012 


0.053 


-0, 


,014 ± 0, 


,012 




1993ag 


0.050 


18.32 


17.83 ± 0, 


,02 


17, 


73 ± 0, 


18 


17.26 ± 0, 


,15 


0.936 ± 0.015 


0, 


,217 ± 0,020 


0.111 


0, 


,120 ± 0, 


,021 


2,3 


1994M 


0.024 


16.34 


16.24 ± 0, 


,03 


16, 


07 ± 0, 


20 


15.84 ± 0, 


16 


0.882 ± 0.015 


0, 


043 ± 0,022 


0.023 


0, 


,063 ± 0, 


,022 




1994S 


0.016 


14.85 


14.78 ± 0, 


,02 


14, 


,83 ± 0, 


,22 


14.86 ± 0, 


,19 


1.033 ± 0.026 


-0, 


,061 ± 0.019 


0.018 


-0, 


,010 ± 0, 


,019 




1995ac 


0.049 


17.23 


17.05 ± 0, 


,01 


17, 


,17 ±0, 


,18 


17.17 ±0, 


,13 


1.083 ± 0.012 


0, 


,026 ± 0.011 


0.042 


-0. 


.005 ± 0, 


,011 




1995bd 


0.016 


17.34 


15.32 ± 0, 


,01 


15, 


,37 ±0, 


,30 






1.039 ± 0.008 


0, 


,735 ± 0.008 


0.490 


0, 


,348 ± 0, 


,009 


1-3 


1996C 


0.030 


16.62 


16.57 ± 0, 


,04 


16, 


,74 ± 0, 


,19 


16.50 ± 0, 


,16 


1.120 ± 0.020 


0, 


,012 ± 0.026 


0.014 


0, 


,051 ± 0, 


,027 




1996ab 


0.125 


19.72 


19.57 ± 0, 


,04 


19, 


,47 ± 0, 


,19 


19.82 ± 0, 


,16 


0.934 ± 0.032 


0, 


,174 ± 0.025 


0.032 


-0, 


,082 ± 0, 


,029 




1996bl 


0.035 


17.08 


16.66 ± 0, 


,01 


16, 


71 ± 0, 


19 


16.55 ± 0, 


,14 


1.031 ± 0.015 


0, 


,093 ± 0.012 


0.099 


0, 


,036 ± 0, 


,012 




1996bo 


0.016 


16.18 


15.85 ± 0, 


,01 


15, 


65 ± 0, 


22 






0.862 ± 0.006 


0, 


,406 ± 0.008 


0.077 


0, 


,383 ± 


,008 


1-3 



a: Supcrnovac through 1993ag arc from HOB, later ones from R,99. 
b: This is the measured peak magnitude of the S-band hghtcurvc. 

c: This includes the Galactic extinction correction and a X-correction: M — B = m^'^°'^ — Kb — Ab, where Kb is the /C-correction and Ab is the Galactic extinction 
correction. The quoted error bar is the uncertainty on the peak magnitude from the lightcurve fit. 

d: This value includes the stretch correction: = m^^"^ — Kb — Ab + (x{s — 1). a is the best-fit value of the stretch/luminosity slope from the fit to the primary 

low-extinction subset (Fit 3 in § 4). The quoted error bar includes all uncertainties for non-extinction corrected fits described in § 2.4. Note that these values are only provided 
for convenience; they were not used directly in any cosmological fits, since the a is also a fit parameter. 

e: Similar to column only with the host-galaxy extinction correction applied. The stretch/luminosity slope used for this value is that from the fit to the primary subset 
(Fit 6 in § 4). The quoted error bar includes all uncertainties for extinction-corrected fits described in § 2.4. A line indicates a supernova which did not appear in the primary 

subset (see § 2.5.) 

/: This value has been iC-corrected and corrected for Galactic extinction. 

g: This is the measured B-V color at the epoch of rest-frame S-band lightcurve maximum. 

h: Schlegel, Finkbeiner, & Davis (1998); this extinction is already included in the quoted values of rriB in column c. 
i: These supernovae are excluded from the indicated subsets; § 2.5. 



tion of time and is described by a spline under 
tension: an "intercept" b{t) and a "slope" m{t). 
At any given date the intrinsic color is 

color{t') = h(i') + m{t') x (l/s^ _ i) (2) 

where t' = t/{s(l + z)), z is the redshift of the 
supernova, and s is the timescale stretch of the 
supernova from a simultaneous fit to the B and V 
lightcurves (matching the procedure used for most 
of the high-redshift supernovae). This arbitrary 
functional form was chosen to match the stretch 
vs. color distribution. 

As the goal was to determine intrinsic colors 
without making any assumptions about redden- 
ing, no host-galaxy extinction corrections were ap- 
plied to the literature data at this stage of the 
analysis. Instead, host-galaxy extinction was han- 
dled by performing a robust blue-side ridge- line fit 
to the supernova color curves, so as to extract the 
unrcddcncd intrinsic color. Individual color points 
that were outliers were prevented from having too 
much weight in the fit with a small added disper- 
sion on each point. The blue ridge- line was se- 
lected by allowing any point more than Ict to the 
red side of the fit model only to contribute to the 

as if it were lo" away. Additionally, those su- 
pernovae which were most reddened were omitted. 
The resulting fit procedure provided B-V ^ V-R^ 
and R-I as a fimction of epoch and stretch; those 
colors were used to correct the template spectrum 
as described above. 

Some of our data extend into the rest-frame V - 
band range of the spectrum. This is obvious for 
supernovae at z > 0.7 where a f7-band template 
is fit to the i?-band data. However, even for su- 
pernovae at z > 0.55, the de-redshifted i?-band 
filter begins to overlap the f7-band range of the 
rest-frame spectrum. Thus, it is also important 
to know the intrinsic U-B color so as to generate 
a proper spectral template. We used data from 
the literature, as given in Table 6. Hero, there is 
an insufficient number of supernova lightcurves to 
reasonably use the sort of ridge-line analysis used 
above to eliminate the effects of host-galaxy ex- 
tinction in determining the intrinsic BVRI colors. 
Instead, for U-B, we perform extinction correc- 
tions using the E(B-V) values from Phillips et al. 
(1999). Based on Table 6, we adopt a, U-B color 
of —0.4 at the epoch of rest-B maximum. This 
value is also consistent with the data shown in 



Jha (2002) for supernovae with timescale stretch of 
s ~ 1, although the data are not determinative. In 
contrast to the other colors, U-B was not consid- 
ered to be a function of stretch. Even though ,Jha 
(2002) does show U-B depending on lightcurve 
stretch, the supernovae in this work that would 
be most affected (those at z > 0.7 where E{B-V) 
is estimated from the rest- frame U -B color) cover 
a small range in stretch; current low-redshift U -B 
data do not show a significant slope within that 
range. See § 5.4 for the effect of systematic error 
in the assumed intrinsic U-B colors. 

Any intrinsic uncertainty in B-V is already sub- 
sumed within the assumed intrinsic dispersion of 
extinction-corrected peak magnitudes (see § 2.4); 
however, we might expect a larger dispersion in 
intrinsic U-B due to e.g., metallicity effects (Hoe- 
flich. Wheeler, & Thielemann 1998; Lentz et al. 
2000). The low-redshift [/-band photometry may 
also have unmodeled scatter e.g., related to the 
lack of extensive UV supernova spectrophotome- 
try for X-corrections. The effect on extinction- 
corrected magnitudes will be further increased by 
the greater effect of dust extinction on the bluer U - 
band light. The scatter of our extinction-corrected 
magnitudes about the best-fit cosmology suggests 
an intrinsic uncertainty in U-B of 0.04 magni- 
tudes. This is also consistent with the U-B data 
of Jha (2002) over the range of timescale stretch of 
our z > 0.7 SNe la, after two extreme color outliers 
from Jha (2002) are removed; there is no evidence 
of such extreme color objects in our dataset. Note 
that this intrinsic U-B dispersion is in addition to 
the intrinsic magnitude dispersion assumed after 
extinction correction. 

The template spectrum which has been con- 
structed may be used to perform color- and K- 
corrections on both the low- and high-redshift su- 
pernovae to be used for cosmology. However, it 
must be further modified to account for the red- 
dening effects of dust extinction in the supernova 
host galaxy, and extinction of the redshifted spec- 
trum due to Galactic dust. To calculate the red- 
dening effects of both Galactic and host-galaxy ex- 
tinction, we used the interstellar extinction law of 
O'Donnell (1994) with the standard value of the 
parameter Ry = 3.1. Color excess {E{B-V)) val- 
ues due to Galactic extinction were obtained from 
Schlegel, Finkbeiner, & Davis (1998). 

The E{B-V) values quoted in Tables 3, 4, and 5 
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Table 6: U-B SN la Colors at Epoch of B-band Maximum 

~SN Raw U-B'' Corrected U-B^ Reference 



1980N 


-0.21 


-0.29 


Hamuy et al. (1991) 


1989B 


0.08 


-0.33 


Wells et al. (1994) 


1990N 


-0.35 


-0.45 


Lira et al. (1998) 


1994D 


-0.50 


-0.52 


Wu, Yan, & Zou (1995) 


1998bu 


-0.23 


-0.51 


Suntzeff et al. (1999) 



a: This is the measured U-B value from the cited paper. 



b: This U-B value is Jf-corrected, and corrected for host-galaxy and 
Galactic extinction. 



are the values necessary to reproduce the observed 
R-I color at the epoch of the maximum of the 
rest-frame B lightcurve. This reproduction was 
performed by modifying the spectral template ex- 
actly as described above, given the intrinsic color 
of the supernova from the fit stretch, the Galac- 
tic extinction, and the host-galaxy E{B-V) pa- 
rameter. The modified spcctnim was integrated 
through the Bessell R- and /-band filters, and 
E{B-V) was varied until the R-I value matched 
the peak color from the lightcurve fit. 

For each supernova, this finally modified spec- 
tral template was integrated through the Bessell 
and WFPC2 filter transmission functions to pro- 
vide color and iiT-corrections. The exact spectral 
template needed for a given data point on a given 
supernova is dependent on parameters of the fit: 
the stretch, the time of each point relative to the 
epoch of rest-B maximum, and the host-galaxy 
E{B-V) (measured as described above). Thus, 
color and /T-corrections were performed iteratively 
with lightcurve fitting in order to generate the fi- 
nal corrections used in the fits described in § 2.2. 
An initial date of maximum, stretch, and host- 
galaxy extinction was assumed in order to gen- 
erate if-corrcctions for the first iteration of the 
fit. The parameters resulting from that fit were 
used to generate new color and JsT-corrections, and 
the whole procedure was rcpeatc;d until the results 
of the fit converged. Generally, the fit converged 
within 2-3 iterations. 

2.4. Cosmological Fit Methodology 

Cosmological fits to the luminosity distance 
modulus equation from the Friedmann- Robertson- 
Walker metric followed the procedure of P99. The 



set of supernova redshifts {z) and i^-corrected 
peak B-magnitudes (ms) were fit to the equation 

mB=M + f>\ogVc{z;VlM,^K)-a{s-l) (3) 

where s is the stretch value for the supernova, 
Vc = H^dL is the "Hubblc-constant-free" lu- 
minosity distance (Perlmutter et al. 1997), and 
M = Mb - 5 logFo -|- 25 is the "Hubble-constant- 
free" _B-band peak absolute magnitude of a s = 1 
SN la with true absolute peak magnitude Mb- 
With this procedure, neither Hq nor Mb need be 
known independently. The peak magnitude of a 
SN la is mildly dependent on the lightcurve decay 
time scale, such that supernovae with a slow de- 
cay (high stretch) tend to be over-luminous, while 
supernovae with a fast decay (low stretch) tend 
to be under-luminous (Phillips et al. 1993); a is a 
slope that parameterizes this relationship. 

There are four parameters in the fit: the mass 
density and cosmological constant f^A, as well 
as the two nuisance parameters, A4 and a. The 
four-dimensional (fijvi, ^^a. a) space is divided 
into a grid, and at each grid point a value is 
calculated by fitting the luminosity distance equa- 
tion to the peak B-band magnitudes and redshifts 
of the supernovae. The range of parameter space 
explored included = [0,3), Cl\ = [—1,3) (for 
fits where host-galaxy extinction corrections are 
not directly applied) or Om = [0,4], Qa = [—1)4) 
(for fits with host-galaxy extinction corrections). 
The two nuisance parameters arc fit in the ranges 
a = [-1,4) and M = [-3.9,3.2). No fur- 
ther constraints are placed on the parameters. 
(These ranges for the four fit parameters contain 
> 99.99% of the probability.) At each point on 
the 4-dimensional grid, a is calculated, and a 
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probability is determined from P oc e"-*^ ' . The 
probability of the whole 4-dimensional grid is nor- 
malized, and then integrated over the two dimen- 
sions corresponding to the "nuisance" parameters. 

For each fit, all peak tub values were cor- 
rected for Galactic extinction using E{B-V) values 
from Schlegel, Finkbeiner, & Davis (1998), using 
the extinction law of O'Donnell (1994) integrated 
through the observed filter.^" For our primary fits, 
the total effective statistical uncertainty on each 
value oirriB included the following contributions: 

• the uncertainty on ms from the lightcurve 
fits; 

• the uncertainty on s, multiplied by a 

• the covariancc between and s; 

• a contribution from the uncertainty in the 
redshift due to peculiar velocity (assumed to 
have a dispersion of 300 km s~^ along the 
line of site); 

• 10% of the Galactic extinction correction; 
and 

• 0.17 magnitudes of intrinsic dispersion 
(H96). 

Fits where host-galaxy extinction corrections are 
explicitly applied use the first five items above 
plus: 

• the uncertainty on E{B-V) multiplied by 
Rb', 

• the covariance between E{B-V) and ms; 

• 0.11 magnitudes of intrinsic dispersion 
(Phillips et al. 1999); and 

• an additional 0.04 magnitudes of intrinsic U- 
B dispersion for z > 0.7. 

Host-galaxy extinction corrections used a value 

Rb s Ab/E{B-V) = 4.1, which results from 
passing a SN la spectrum through the standard 
O'Donnell (1994) extinction law. Except where 
explicitly noted below, the E{B-V) uncertainties 
are not reduced by any prior assumptions on the 
intrinsic color excess distribution. Although there 
is almost certainly some intrinsic dispersion cither 
in Rb, or in the true B-V color of a SN la (Nobili 
et al. 2003), we do not explicitly include such a 

'This supersedes P99, where an incorrect dependence on z 
of the effective Rfj for Galactic extinction was applied. The 
corrected procedure decreases the flat- universe value of Um 
by 0.03. 



term. The effect of sticli a dispersion is included, 
in principle, in the 0.11 magnitudes of intrinsic 
magnitude dispersion which Phillips et al. (1999) 
found after applying extinction corrections. 

As discussed in § 2.3, the intrinsic U-B disper- 
sion is likely to be greater than the intrinsic B-V 
dispersion. For those supernovae most affected by 
this (i.e. those at 2; > 0.7), we included an ad- 
ditional uncertainty corresponding to 0.04 magni- 
tudes of intrinsic U-B dispersion, converted into 
a magnitude error using the O'Donnell extinction 
law. 

This set of statistical uncertainties is slightly 

different from that used in P99. For these fits, the 
test value of a was used to propagate the stretch 
errors into the corrected S-band magnitude er- 
rors; in contrast, P99 used a single value of a for 
purposes of error propagation. 

2.5. Supernova Subsets 

In P99, separate analyses were performed and 
compared for the supernova sample before and af- 
ter removing supernovae with less secure identifi- 
cation as Type la. The results were shown to be 
consistent, providing a cross-check of the cosmo- 
logical conclusions. For the analyses of this paper, 
adding and comparing eleven very-wcU-measured 
SNe la, we only consider from P99 the more se- 
curely spectrally identified SNe la with reason- 
able color measurements (i.e. (Jr-i < 0.25); those 
supernovae are listed in Table 4. Following P99, 
we omit one supernova which is an outlier in the 
stretch distribution, with s < 0.7 (SN 1992br), and 
one SN which is a > 6cr outlier from the best- 
fit cosmology (SN 19970). We also omit those 
supernovae which are most seriously reddened, 
with E{B-V) > 0.25 and > 3a above zero; host- 
galaxy extinction corrections have been found in 
studies of low-redshift supernovae to overcorrect 
these reddest objects (Phillips et al. 1999). This 
cut removes two SNe at low redshift (SNe 1995bd 
and 1996bo), one from P99 (SN1996cn), and one 
of the eleven HST supernovae from this paper 
(SN1998aw). The resulting "full primary subset" 
of SNe la is identified as Subset 1 in the tables. 

For the analyses of a "low-extinction primary 
subset," Subset 2, we further cull out four super- 
novae with host-galaxy E{B-V)> 0.1 and > 2a 
above zero, including two of the HST supernovae 
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from this paper (SNc 1992ag, 1993ag, 1998as, and 
1998ax). The low-extinction primary subset in- 
cludes eight of the eleven new HST supernovae 
presented in this paper. 

Subset 3, the "low-extinction strict la subset," 
makes an even more stringent cut on spectral con- 
firmation, including only those supernovae whose 
confirmations as Type la SNe are unquestionable. 
This subset is used in § 5.2 to estimate any pos- 
sible systematic bias resulting from type contam- 
ination. An additional six supernovae, including 
two of the HST supernovae from this paper, are 
omitted from Subset 3 beyond those omitted from 
Subset 2; these are SNel995as, 1996cf, 1996cg, 
1996cm, 1998ay, and 1998be. 

3. Colors and Extinction 

In this section, we discuss the limits on host- 
galaxy extinction we can set based on the mea- 
sured colors of our supernovae. For the primary 
fit of our P99 analysis, extinction was estimated 
by comparing the mean host-galaxy E{B-V) val- 
ues from the low- and high-redshift samples. Al- 
though the uncertainties on individual E{B-V) 
values for high-redshift supernovae were large, the 
uncertainty on the mean of the distribution was 
only 0.02 magnitudes. P99 showed that there was 
no significant difference in the mean host-galaxy 
reddening between the low and high-redshift sam- 
ples of supernovae of the primary analysis (Fit 
C). This tightly constrained the systematic un- 
certainty on the cosmological results due to dif- 
ferences in extinction. The models of Hatano, 
Branch, & Deaton (1998) suggest that most SNe la 
should be found with little or no host galaxy ex- 
tinction. By making a cut to include only those 
objects which have small E{B-V) values (and 
then verifying the consistency of low- and high- 
redshift mean reddening), we are creating a sub- 
sample likely to have quite low extinction. The 
strength of this method is that it does not de- 
pend on the exact shape of the intrinsic extinction 
distribution, but only requires that most super- 
novae show low extinction. Figure 3 (discussed 
below) demonstrates that most supernovae indeed 
have low-extinction, as expected from the Hatano, 
Branch, & Deaton (1998) models. Monte Carlo 
simulations of our data using the Hatano, Branch. 
& Deaton (1998) extinction distribution function 



and our low-extinction E{B-V) cuts confirm the 
robustness of this approach, and further, demon- 
strate that similarly low extinction subsamples are 
obtained for both low- and high-redshift datasets 
despite the larger color uncertainties for some of 
the P99 supernovae. 

Riess et al. (1998) used the work of Hatano, 
Branch, & Deaton (1998) difi^erently, by apply- 
ing a one-sided Bayesian prior to their measured 
E{B-V) values and uncertainties. A prior formed 
from the Hatano, Branch, & Deaton (1998) extinc- 
tion distribution function would have zero proba- 
bility for negative values of E{B-V), a peak at 
E{B-V) ~ with roughly 50% of the probabil- 
ity, and an exponential tail to higher extinctions. 
As discussed in P99 (see the "Fit E" discussion, 
where P99 apply the same method), when uncer- 
tainties on high- and low-redshift supernova col- 
ors differ, use of an asymmetric prior may intro- 
duce bias into the cosmological results, depending 
on the details of the prior. While a prior with a 
tight enough peak at low extinction values intro- 
duces little bias (especially when low- and high- 
redshift supernovae have comparable uncertain- 
ties), it does reduce the apparent E{B-V) error 
bars on all but the most reddened supernovae. As 
we will show in Figure 10 (§ 4.1) the use of this 
prior almost completely eliminates the contribu- 
tion of color uncertainties to the size of the cosmo- 
logical confidence regions, meaning that an extinc- 
tion correction using a sharp enough prior is much 
more akin to simply selecting a low-extinction sub- 
set than to performing an assumption-free extinc- 
tion correction using the E{B-V) measurement 
uncertainties. 

The high precision measurements of the R-I 
color afforded by the WFPC2 lightcurves for the 
new siipernovae in this work allow a direct estima- 
tion of the host-galaxy E{B-V) color excess with- 
out any need to resort to any prior assumptions 
concerning the intrinsic extinction distribution. 

Figure 3 shows histograms of the host-galaxy 
E{B-V) values from different samples of the su- 
pernovae used in this paper. For the bottom 
two panels, a line is over-plotted that treats the 
H96 low-extinction subset's E{B-V) values as a 
parent distribution, and shows the expected dis- 
tribution for the other samples given their mea- 
surement uncertainties. The low-extinction sub- 
set of each sample (the grey histogram) has a 
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Fig. 3. — Histograms of E{B-V) for the 
four samples of supernovae used in this 
paper. The filled grey histogram repre- 
sents jtist the low-extinction subset (Stib- 
set 2). The open boxes on top of that 
represent supernovae which are in the pri- 
mary subset (Subset 1) but excluded from 
the low-extinction subset. Finally, the 
dotted histogram represents those super- 
novae which are in the ftill sample but 
omitted from the primary subset. The 
solid lines drawn over the bottom two 
panels is a simulation of the distribution 
expected if the low-extinction subset of 
the H96 sample represented the true dis- 
tribution of SN colors, given the error bars 
of the low-extinction subset of each high- 
redshift sample. 
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Table 7: Mean E{B-V) Values 



Sample 


Complete 


Low-extinction 




Set 


Primary Subset 






SNe" 


Low z 


+0.095 ± 0.003 


-0.001 + 0.003 


P99 


+0.018 ±0.024 


-0.004 + 0.025 


HST 


+0.090 + 0.012 


+0.012 + 0.015 



a: SNe omitted from our low-extinction pri- 
mary subset, Subset 2, (§ 2.5) have been omit- 
ted from these means. This excludes outliers, as 
well as supernovae with both E{B-V) > 0.1 and 
E{B-V) > 2a above zero. 



color excess distribution which is consistent with 
that of the low-extinction subset of H96. Table 7 
lists the variance-weighted mean E{B-V) values 
for the low-redshift supernovae and for each sam- 
ple of high-redshift supernovae. Although varying 
amounts of extinction arc detectable in the mean 
colors of each full sample, the supernovae in the 
low-extinction primary subset (§ 2.5) of each sam- 
ple are consistent with E{B-V) = 0. This subset 
is consistent with the models of Hatano, Branch, 
& Deaton (1998), discussed above, in which most 
SNe la arc observed in regions of very low ex- 
tinction. We will consider cosmological fits both 
to this low-extinction subset and to the primary 
subset with host-galaxy reddening corrections ap- 
plied. 

Figure 4 shows E{B-V) vs. redshift for the 
eleven supernovae of this paper. Three of the 
lowest redshift SNe are likely to be significantly 
reddened: SN 1998as at z = 0.36, SN 1998aw at 
z = 0.44, and SN 1998ax at z = 0.50. This higher 
incidence of extincted SNe at the low-redshift end 
of our sample is consistent with expectations for 
a flux-limited survey, where extincted supernovae 
will be preferentially detected at lower redshifts. 
Indeed, the distribution of E{B-V) values versus 
redshift shown in Figure 4 is consistent with the 
results of a Monte Carlo simulation similar to that 
of Hatano, Branch, & Deaton (1998), but includ- 
ing the effects of the survey flux limit. Several au- 
thors (including Leibundgut (2001) and Falco et 
al. (1999)) have suggested that there is evidence 
from the E{B-V) values in Riess et al. (1998) 
that high-redshift supernovae are bluer statisti- 
cally than their low-redshift counterparts. Our 
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Fig. 4. — A plot of E{B-V) as a function of redshift 

for the 11 HST-obscrvcd supernovae of this paper 
shows that the blue edge of the distribution shows 
no significant evolution with redshift. (The larger 
dispersion at lower redshifts is expected for a flux- 
limited sample.) Error bars include only measure- 
ment errors, and no assumed intrinsic color disper- 
sion. Filled circles are those supernovae in the low- 
extinction subset (Subset 2). 



data show no such effect (nor did our P99 SNe). 

The mean host-galaxy color excess calculated 
for the highest redshift supernovae is critically de- 
pendent on the assumed intrinsic U-B color (see 
§ 2.3). An offset in this assumed U-B will af- 
fect the high-redshift supernovae much more than 
the low-redshift supernovae (whose measiirements 
are primarily of the rest frame B- and F-band 
lightcurves). The /f-corrected, rest-frame B-band 
magnitudes arc also dependent on the assumed 
supernova colors that went into deriving the K- 
corrections. If the assumed U-B color is too red, 
it will affect the cross-filter X-correction applied 
to i?-band data at z > 0.5, thereby changing de- 
rived rest frame colors. In § 5, we consider the 
effect of changing the reference U-B color. 

4. Cosmological Results 
4.1. and ilj^ 

Figures 5 through 7 show Hubble Diagrams 
which effective B-band peak magnitudes and red- 
shifts for the new supernovae of this paper; these 
magnitudes have been K- and stretch corrected, 
and have been corrected for Galactic extinction. 
Figure 5 shows all of the data in the low-extinction 
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Fig. 5. — Hubble diagram of effective K- and stretcli-corrected niB vs. redshift for the supernovae 
in tlie primary low-extinction subset. Filled circles represent the HST supernovae of this paper. 
Inner error bars show just the measurement uncertainties; outer error bars include 0.17 magnitudes 
of intrinsic dispersion. The solid line is the best-fit flat-universe cosmology from the low-extinction 
subset; the dashed and dotted lines represent the indicated cosmologies. 
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Fig. 6. — Upper panel: Averaged Hubble diagram with a linear redshift scale 
for all supernovae from our low-extinction subsample. Here supernovae within 
Az < 0.01 of each other have been combined using a weighted average in order 
to more clearly show the quality and behavior of the dataset. (Note that these 
averaged points are for display only, and have not been used for any quantitative 
analyses.) The solid curve overlaid on the data represents our best-fit flat-universe 
model, {VLm.^a) = (0.25,0.75) (Fit 3 of Table 8). Two other cosmological mod- 
els are shown for comparison: (f^Mj^^A) = (0.25,0) and (f^Mi^^A) = (1,0). Lower 
panel: Residuals of the averaged data relative to an empty universe, illustrating the 
strength with which dark energy has been detected. Also shown are the suite of 
models from the upper panel, including a solid curve for our best-fit flat-universe 
model. 23 
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Fig. 7. — Hubble diagram of effective K- 

and strctch-corrcctcd vs. rcdshift for 
the 11 supernovae observed with WFPC2 
and reported in this paper. Circles rep- 
resent supernovae in the primary subset 
(Subset 1); the one point plotted as a cross 
(the very reddened supernova SN 1998aw) 
is omitted from that subset. Open cir- 
cles represent reddened supernovae omit- 
ted from the low-extinction primary sub- 
set (Subset 2), while filled circles arc in 
both Subsets 1 and 2. Upper plot: 
no host-galaxy E{B-V) extinction cor- 
rections have been applied. Inner error 
bars only include the measurement er- 
ror. Outer error bars include 0.17 mag- 
nitudes of intrinsic dispersion. Lower 
plot: extinction corrections have been ap- 
plied using the standard interstellar ex- 
tinction law. Error bars have been in- 
creased by the uncertainty in this ex- 
tinction correction. Again, inner error 
bars represent only measurement \mcer- 
tainties, while outer error bars include 
0.11 magnitudes of intrinsic dispersion. 
Lines are for three example cosmologies 
with the indicated values of Qm and f2A; 
the solid line is the best-fit flat-universe 
cosmology to our full primary subset with 
extinction corrections applied. 
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Fig. 8.— 68%, 90%, 95%, and 99% confidence regions for flu and 
from this paper's primary analysis, the fit to the low-extinction primary 
subset (Fit 3). 
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subset of supcrnovac. For the sake of clarity, Fig- 
ure 6 shows the same subset, but for this figure su- 
pernovae with redshifts within 0.01 of each other 
have been combined in a variance-weighted aver- 
age. The lower panel of Figure 6 shows the residu- 
als from an empty universe (flu = 0, Oa = 0), il- 
lustrating the strength with which dark energy has 
been detected. In both Figures 5 and 6, the solid 
line represents the flat-universe cosmology result- 
ing from our fits to the low-extinction subset. Fig- 
ure 7 shows just the eleven HST supernovae from 
this paper. In the upper panel of this latter figure, 
the stretch- and X-corrected effective uib values 
and uncertainties are plotted. In the lower panel, 
eff'ective tub values have also been corrected for 
host-galaxy extinction based on measured E{B-V) 
values. The solid line in this figure represents the 
best-fit flat-universe cosmology to the full primary 
subset with extinction corrections applied. 

Table 8 lists results from fits to both of our 
primary subsets of supernovae. Supernovae from 
both the H96 and R99 low-redshift samples were 
included in all fits. The first three lines show fits 
to the low-extinction primary subset. So that the 
new sample of high-redshift supernovae may be 
compared to those from P99, each high-redshift 
sample was fit separately (Fits 1 and 2). Fit 3 
combines all of the current SCP high-redshift su- 
pernovae from the low-extinction subsets, and rep- 
resents the primary result on Om and Qa for this 
paper; Figure 8 shows the confidence regions for 
f^M vs. Oa from this fit. Figure 9 shows the com- 
parison of the confidence regions when each high- 
redshift sample is treated separately. Note that 
Fit 2 provides comparable and consistent measure- 
ments of CIm and to Fit 1. Additionally, the 
sizes of the confidence regions from the 8 HST SNc 
in Fit 2 is similar to those in Fit 1 , which includes 
25 high-redshift supernovae from P99. 

Fits 4-6 in Table 8 show the results for the 
primary subset when host-galaxy extinction cor- 
rections have been applied. Figure 10 compares 
these results to those of the primary low-extinction 
fit. The primary fits of Figure 9 are reproduced 
in the top row of Figure 10. The second row 
has host-galaxy extinction corrections applied us- 
ing the one-sided prior used by Fit E of P99 and 
Riess et al. (1998) discussed in § 3. The third 
row has full extinction corrections applied with- 
out any prior assumptions on the intrinsic E[B-V) 
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Fig. 9.— Contours indicate 68% and 90% 
confidence regions for fits to supernovae 
from the low-extinction primary subset, 
including just the high-redshift SNe from 
P99 (dotted lines), just the new HST 
high-redshift SNe (solid lines), and all 
SCP high-redshift SNe (filled contours). 
The low-redshift SNe from the primary 
subset are included in all fits. The new, 
independent sample of high-redshift su- 
pernovae provide measurements of Q.m 
and r^A consistent with those from the 
P99 sample. 



distribution. Three conclusions are apparent from 
this plot. First, using a strongly peaked prior on 
extinction prevents the E{B-V) error bars from 
being fully propagated into the cosmological con- 
fidence regions, and hence apparently tightens the 
constraints. However, for a peaked prior, this is 
very similar to assuming no extinction and not 
performing an extinction correction (but without 
testing the assumption), while for a wider prior 
there is a danger of introducing bias. Second, the 
current set of supernovae provide much smaller 
confidence regions on the f^A versus plane than 
do the SNe la from previous high-redshift sam- 
ples when unbiased extinction corrections are ap- 
plied. Whereas Figure 9 shows that the current 
set of supernovae give comparable measurements 
of Q.M and Q.^ when the low-extinction subsample 
is used with no host-galaxy extinction corrections. 
Figure 10 shows that the much higher precision 



Table 8: Cosmological fits 



Fit High-Redshift SNe 


NsNe Min. 


Hm for 


Q\ for 


P(Qa > 0) 


M 


a 


# Included in Fit" 




Flat'' 


Flat'' 








Fits to the Low-Extinction Primary Subset 












1 SNe from P99 


46 52 


n 25+° °^ 


"•'^-0.07 


0.9995 


-3.49 ± 0.05 


1.58 ±0.31 


2 New HST SNe 
from this paper 


29 30 


U-^O_0.08 


„ 7^+0.08 


0.9947 


-3.47 ±0.05 


1.06 ± 0.37 


3 All SCP SNe 


54 60 


'-'•■'u-o.oe 


'"-0.07 


0.9997 


-3.48 ±0.05 


1.47 ± 0.29 


Fits to Full Primary Subset, with Extinction Correction 










4 SNe from P99 


48 56 




"■'^-0.18 


0.9967 


-3.55 ±0.05 


1.30±0.30 


5 New HST SNe 
from this paper 


33 39 




n 70+0.10 


0.9953 


-3.54 ±0.05 


1.29±0.28 


6 All SCP SNe 


58 65 




'^-0.11 


0.9974 


-3.53 ± 0.05 


1.18±0.30 


a: All fits include the low-redshift SNe from H96 and R99. 


See § 2.5 for the definitions of the supernova 


subsets. 


b: This is the intersection of the fit probability distribution with the line CIm + = 


1. 





color measurements from the WFPC2 data allow 
us directly to set much better limits on the efii'ects 
of host-galaxy extinction on the cosmological re- 
sults. Finally, the cosmology which results from 
the extinction-corrected fits is consistent with the 
fits to our low-extinction primary subset. Con- 
trary to the assertion of Rowan-Robinson (2002), 
even when host-galaxy extinction is directly and 
fully accounted for, dark energy is required with 
P(Oa > 0) > 0.99. 

4.2. Combined High-Redshift Supernova 
Measurements 

Figure 11 shows measurements of SIm and Qa 
which combine the high-rcdshift supernova data 
of Riess et al. (1998) together with the SCP data 
presented in this paper and in P99. The con- 
tours show confidence intervals from the 54 su- 
pernovae of the low-extinction primary Subset 2 
(used in Fit 3 of Table 8), plus the nine well- 
observed confirmed Type la supcrnovac from Riess 
et al. (1998) (using the lightcurve parameters re- 
sulting from their template- fitting analysis); fol- 
lowing the criteria of Subset 2, SN 1997ck from 
that paper has been omitted, as that supernova 



was not confirmed spectrally. We also omit from 
Riess et al. (1998) the supernovae they measured 
using the "snapshot" method (due to the very 
sparsely sampled lightcurve), and two SCP su- 
pcrnovac that Riess et al. (1998) used from the 
P99 data set which are redundant with our sam- 
ple. This fit has a minimum of 65 with 63 
supernovae. Under the assumption of a flat uni- 
verse, it yields a measurement of the mass density 
of Om = 0.26to;o6' equivalently a cosmological 
constant of = 0.74^°:°?. Recent ground-based 
data on eight new high-redshift supernovae from 
Tonry et al. (2003) (not included in this fit) are 
consistent with these results. Note that in this 
fit, the nine supernovae from Riess et al. (1998) 
were not treated in exactly the same manner as 
the others. The details of the template fitting will 
naturally have been different, which can introduce 
small differences (see § 5.1). More importantly, 
the iiT-corrections applied by Riess et al. (1998) to 
derive distance moduli were almost certainly dif- 
ferent from those used in this paper. 



27 



3r ' ' ' M " ■ I ' " ' 'J L' ' ' ' r ' ' ' I • ' I '11 " ' I ' ' ' ' I ' ' ' ' 




Fig. 10. 68.3%, 95.4%. and 99.7% confidence regions for Qm and Ha using different data subsets and methods for treating 

host-galaxy extinction corrections. The top row represents our fits to the low-extinction primary subset, where significantly reddened 
supernovae have been omitted and host-galaxy extinction corrections are not applied. The second row shows fits where extinction 
corrections have been applied using a one-sided extinction prior. These fits are sensitive to the choice of prior, and can cither yield 
results equivalent to analyses assuming low extinction (but without testing the assumption), or yield biased results (see text). Note that 
the published contours from Riess et al. (1998; their Fig. 6, solid contours) presented results from fits that included nine well-observed 
supernovae (that are comparable to the primary subsets used in the other panels), but also four supernovae with very sparsely sampled 
lightcurves, one supernova at z = 0.97 without a spectral confirmation, as well as two supernovae from the P99 set. The third row 
shows fits with unbiased extinction corrections applied to our primary subset. The HST SNe presented in this paper show a marked 
improvement in the precision of the color measurements, and hence in the precision of the Qm and Qa measurements when a full 
extinction correction is applied. With full and unbiased extinction corrections, dark energy is still required with P{Qa > 0) = 0.99. 
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Fig. 11.— 68%, 90%, 95%, 
and 99% confidence regions for 
flu and Q\, combining the high- 
redshift data of the SCP (this 
paper and P99) and Riess et al. 
(1998). The fit includes Subset 2 
supernovae from the SCP plus 
the nine well-observed confirmed 
SNe la from Riess et al. (1998). 



4.3. Dark Energy Equation of State 

The fits of the previous section used a tradi- 
tional constrained cosmology where CIm is the en- 
ergy density of non-relativistic matter (i.e. pres- 
sure p = 0), and CIa is the energy density in a cos- 
mological constant (i.e. pressure p = —p, where 
p is the energy density). In Einstein's field equa- 
tions, the gravitational effect enters in terms of 
p + 3p. If ID = p/ p is the equation of state pa- 
rameter, then for matter w = 0, while for vac- 
uum energy (i.e. a cosmological constant) w = —1. 
In fact, it is possible to achieve an accelerating 
Universe so long as there is a component with 
w <~ — 1/2. (If there were no contribution from 
f^M, only w < —1/3 dark energy is necessary for 
acceleration; however, for plausible mass densities 
^ 0.2, the dark energy must have a more neg- 
ative value of u>.) The Hubble diagram for high- 
redshift supernovae provides a measurement of w 
(P99, Garnavich et al. 1998b). Panels (a) and (b) 
of Figure 12 show the joint confidence regions for 
f^M versus w from the SCP supernovae, including 



the new HST supernovae, under the assumptions 
that w is constant with time, and that the Uni- 
verse is flat, i.e. + = 1 (where fix is the 
energy density in the component with equation of 
state w, in units of the critical density). The su- 
pernova alone data set a 99% confidence limit of 
w < —0.64 for any positive value of flu, without 
any prior assumptions on w. 

A fit with extinction corrections applied to the 
full primary subset (Fit 6, shown in Figure 12b) 
gives a 99% confidence limit of w < —1.00. How- 
ever, this latter limit should be approached with 
caution, because w is not well bounded from be- 
low with the supernova data alone. Although Fig- 
ure 12 only shows confidence intervals down to 
w = —2, the 68% confidence interval from Fit 3 
extends to w < —4, and the 99% confidence inter- 
val extends tow < —10; these confidence intervals 
extend to even further negative w in Fit 6. The 
weight of probability at very low (and probably 
implausible) w pulls the 68% confidence interval 
in Fit 6 (Figure 12b) downward. A fit which used 
a prior to restrict w to more reasonable values (say 
w > —2) would show similar outer confidence in- 
tervals, but a 68% confidence interval more similar 
to that of the low-extinction subset in Figure 12a. 

Other methods provide measurements of flu 
and w which are complementary to the supernova 
results. Two of these measurements are plotted 
in the middle row of Figure 12, compared with 
the supernova measurements (in dotted contours) . 
In filled contours are results from the redshift- 
distortion parameter and bias-factor measurement 
of the 2dF Galaxy Rcdshift Survey (2dFGRS) 
(Hawkins et al. 2002; Verde et al. 2002). These 
provide a measurement of the growth parameter, 
/ = 0.51 ± 0.11, at the survey rcdshift z = 0.15. 
We have used the method of Linder & Jenkins 
(2003) to directly solve for /(Om, w, z) rather than 
convert / to Hm, as the conversion formula given 
in Hawkins et al. (2002) is valid only for w = —1. 
Comparison of the 2dFGRS value of / with the 
calculated values of f{flM,w,z) yields the joint 
confidence region for 17m and w.^^ 

In solid lines in panels (c) and (d) of Fig- 
ure 12 are contours representing confidence re- 



Note that we have not used the independent 2dFGRS power 

spectrum constraint on because it has not yet been 

generalized for diflferent values of w. 
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Fig. 12. — Joint measurements of flu and w assuming f^M + = 1 and that w is not time-varying. 
Confidence regions plotted are 68%, 90%, 95%, and 99%. The left column (panels a, c, and e) shows fits to 
the low-extinction primary subset; the right column (panels b, d, and f) shows fits to the primary subset 
with unbiased individual host-galaxy extinction corrections applied to each supernova. The upper panels (a 
and b) show the confidence intervals from the SCP supernovae alone. The middle panels (c and d) overlay 
this (dotted lines) with measurements from 2dFGRS (filled contours) (Hawkins et al. 2002) and combined 
CMB measurements (solid contours) (Bennett et al. 2003; Spergel et al. 2003). The bottom panels (e and f) 
combine the three confidence regions to provide a combined measurement of flu and w. 
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gions basc'd on the distance to the surface of last 
scattering &t z = 1089 from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) and other 
CMB measurements (Bennett et al. 2003; Spcrgcl 
et al. 2003). For a given and w, this reduced 
distance to the surface of last scattering, 7, is given 
by: 

/•1089 

Jo 

(4) 

The plotted CMB constraints come from the 

"WMAPext" sample, which includes other CMB 
experiments in addition to WMAP. They yield a 
measurement of Iq = 1.76 ± 0.058, corresponding 
to r^M = 0.29 at w = —1. Confidence intervals 
are generated by calculating a probability using 
= [{I — Io)/cTig]^, where I is calculated for each 

f^M, W. 

As both of these measurements show mild cor- 
relations between flu and it; in a different sense 
from that of the supernova measurement, the com- 
bined measurements provide much tighter over- 
all constraints on both parameters. The con- 
fidence regions which combine these three mea- 
surements are shown in panels (e) and (f) of 
Figure 12. When the resulting probability dis- 
tribution is marginalized over flu, we obtain a 
measurement oi w = —1.05fo[lo (for the low- 
extinction subset), or w; = — I.O2+Q24 (for the full 
primary subset with host-galaxy extinc;tion cor- 
rections applied). When the probability distri- 
bution is marginalized over w, we obtain a flat- 
universe measurement of flu = 0.27^q'q5 (for the 
low-extinction subset), or Qm — 0.28lo;o5 (for 
the primary subset with host-galaxy extinction 
corrections applied). The 95% confidence limits 
on w when our data are combined with CMB 
and 2dFGRS are -1.6K w < -0.78 for the low- 
extinction primary subset, or —1.67 < m; < —0.62 
for the full extinction-corrected primary subset. If 
we add an additional prior that w > —1, we ob- 
tain a 95% upper confidence Hmit of w < —0.78 for 
the low-extinction primary subset, ot w < —0.67 
for the extinction-corrected full primary subset. 
These values may be compared with the limit in 
Spergel et al. (2003) which combines the CMB, 
2dFGRS power spectrum, and HST key project 
Hq measurements to yield a 95% upper limit of 
w < —0.78 assuming w > —1. Although both our 



measurement and that of Spcrgcl et al. (2003) in- 
clude CMB data, they are complementary in that 
our limit does not include the Hq prior, nor does it 
include any of the same external constraints, such 
as those from large scale structure. 

These combined measurements remain consis- 
tent with a low density universe dominated by 
vacuum energy (constant w = —1), but are also 
consistent with a wide range of other both time- 
varying- w and constant- w dark energy models. 

5. Systematic Errors 

The effect of most systematic errors in the 
vs. Qa plane is asymmetric in a manner similar to 
the asymmetry of our statistical errors. For the ef- 
fects listed below, a systematic difference will tend 
to move the confidence ellipses primarily along 
their major axis. In other words, these systematic 
effects produce a larger uncertainty in I^m + 
than in Om — ^^A (or, equivalently, in a measure- 
ment of Qm or Cl\ alone under the assumption of a 
flat universe). This means that systematic effects 
do not currently hamper the cosmological mea- 
surements from supernovae where they have the 
greatest weight relative to other techniques, nor 
do they significantly diminish the direct evidence 
from supernovae for the presence of dark energy. 
However, they do limit the ability of supernovae to 
measure the spatial curvature ("geometry") of the 
Universe. (Note that the semi-major axis is not 
precisely in the direction of ^lM + ^A, nor is the 
semi- minor axis precisely aligned with Qm — ^A, 
but since these are useful constraints we will quan- 
tify the systematic uncertainties along these two 
directions.) Figure 13 shows the effects of some of 
the systematics discussed in the following subsec- 
tions. 

Systematic effects on flat-universe measure- 
ments of w arc smaller than the current statis- 
tical uncertainties. The right column of Figure 13 
shows the effect of the systematics on the SIm 
versus w confidence regions derived from our su- 
pernova data alone. To quantify the effect of 
identified systematics in the following subsections, 
we determine the shift in the maximum-likelihood 
value of w when the supernova data is combined 
with the CIm versus w confidence regions from 2dF- 
GRS and the CMB (See § 4.3.) 
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Fig. 13.— Simulated effects of iden- 
tified systematic errors on the cosmologi- 
eal parameters, estimated by applying the 
systematic effect to the supernova param- 
eters used in the cosmological fits. The 
left column shows fits to SIm and 0,a, and 
the right column to Qm and the dark en- 
ergy equation of state parameter w. Rows 
{a)-(c) show our primary fit (Fit 3) in 
filled contours, (a) The dotted contours 
show the resiilts of a fit to Subset 3, 
only those supernovae with the most se- 
cure spectral identifications as Type la 
SNe. (b) The dotted contours show a fit 
to Subset 1 where the supernova magni- 
tudes have been dimmed to correct for 
Malmquist bias. (c) The dotted con- 
tours show a fit to Subset 2, where K- 
corrections have been applied using a tem- 
plate spectrum with an intrinsic value of 
U-B— — 0.5 at the epoch of B-maximum. 
(d) The filled contours is Fit 6, the fit to 
the full primary subset with host-galaxy 
extinction corrections applied; the dot- 
ted contours show a fit to the same Sub- 
set, but using a template spectrum with 
an intrinsic value of U-B— — O.C} for esti- 
mating both /^-corrections and color ex- 
cesses, (c) The dotted contours apply ex- 
tinction corrections to Subset 1 using a 
value of Rb = 3.5 rather than the stan- 
dard Rb =4.1 which was used for Fit 6 
(filled contours). 
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5.1. Fit Method, Subset Selection, and 

Choice of a 

There are multiple reasonable choices for 
lightcurve fitting methods which yield slightly dif- 
ferent results for the lightcurve parameters. For 
the supernovae in P99, the i?-band data on high- 
redshift supernovae provided much stronger limits 
on the stretch (the shape of the lightcurve) than 
did more sparse /-band lightcurves. For consis- 
tency, in P99 the stretch values for the low-redshift 
supernovae wore therefore measured using only the 
i?-band lightcurves. 

In this paper, there are high-quality photomet- 
ric measurements from WFPC2 in both R and / 
bands. Thus, data in both colors contribute signif- 
icantly to the constraints on stretch. Additionally, 
photometry is extracted from HST and ground- 
based images in very different apertures, mean- 
ing that different amounts of host galaxy light will 
be included; this background must be subtracted 
from each before the two are combined. As such, 
it is more appropriate to fit these supernovae with 
fixed rather than floating lightcurve zero offsets. 
As this is the most appropriate fit method for 
the HST supernovae, the low-redshift supernovae 
should be treated consistently. These procedures 
which are most appropriate for the HST super- 
novae were used for all new fits performed in this 
paper, and listed in Tables 3 through 5. 

To estimate the size of the effect due to these 
differences in fitting method, cosmological confi- 
dence intervals were generated from the "Case C" 
subset of P99 using the new fits presented in this 
paper and compared to the results quoted in P99 
and other variations on the fitting method. Dif- 
ferences in the fit method can change the flat- 
universe value of f^M by ~0.03, and the value of 
f^M + by up to ~0.8. (This is still much less 
than the major-axis extent of the statistical confi- 
dence ellipse in this direction.) We use these val- 
ues as the "fit-method" systematic uncertainties. 
We similarly performed joint fits to U,yi,w in the 
flat-universe, constant-w case to the supernovae 
from P99 with different lightcurve fit methodolo- 
gies, and from these fits we adopt a fit-method 
systematic uncertainty of 0.02 on constant w (once 
combined with measurements from 2dFGRS and 
the CMB). 

We have also performed a fit without any 



stretch correction at all, i.e. using fixed a Al- 
though the quality of the fit is worse (x^ = 82 with 
54 supernovae, in comparison to = 60 from 
Fit 3), it yields consistent cosmological results, 
with shifts (AQ^* < 0.01) much smaller than the 
already- adopted "fit method" systematic. We 
have likewise performed a fit to the complete set of 
supernovae (including all from P99 with measured 
colors). The fit cosmological values are similarly 
consistent with the primary low-extinction fit. We 
therefore conclude that the effects of these choices 
are subsumed in the "fit method" systematic. 

5.2. Non-Type la Supernova Contamina- 
tion 

All subsets of supernovae used for cosmologi- 
cal fits in this paper omit supernovae for which 
there is not a spectral confirmation of the super- 
nova type. Nonetheless, it is possible in some cases 
where that confirmation is weak that we may have 
contamination from non-Type la supernovae. To 
estimate such an effect, we performed fits using 
only those supernovae which have a firm identi- 
fication as Type la; this is the "strict-la subset" 
from § 2.5. The comparison between our primary 
fit (Fit 3) and this fit with a more stringent type 
cut is shown in row (a) of Figure 13. This fit has a 
value of Om in a flat universe which is 0.03 higher 
than that of Fit 3. The value of Qm + is 0.48 
lower than that of Fit 3. We adopt these values as 
our "type contamination" systematic error. 

The size of this systematic for w is shown in the 
right panel of Figure 13a. Combined with CMB 
and 2dFGRS measurements, the best-fit value of 
w is larger by 0.07; we adopt this as our type con- 
tamination systematic error on w. 

5.3. Malmquist Bias 

As most of our supernovae are from flux- 
limited samples, they will suffer Malmquist bias 
(Malmquist 1924, 1936). This effect was discussed 
extensively in P99, and here we update that dis- 
cussion to include our new HST SNe la. For the 
measurement of the cosmological parameters, it 
is the difference between the Malmqiiist bias of 
the low-redshift and high-redshift samples which 
matters. In particular, the apparent probability 
of > is enhanced only if the low-redshift 
supernovae suffer more Malmquist bias than the 
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high-rcdshift supcrnovac, as this makes the high- 
redshift SNe la seem fainter. 

The P99 high-redshift dataset was estimated to 
have httle Malmquist bias (0.01 mag) because the 
SN discovery magnitudes were decorrelated with 
the measured peak magnitudes. However, for the 
new HST sample, nine of the eleven SNe la (se- 
lected from larger samples of supernovae found in 
the searches) were found almost exactly at max- 
imum Ught. This may reflect a spectroscopic 
flux limit superimposed on the original search flux 
limit since only spectroscopically confirmed SNe la 
were considered, and of those, generally the higher 
redshift SNe la from a given search were chosen 
for HST for follow-up. In particular, the SNe la 
selected for follow-up from the fall 1997 search 
were all found at maximum light, while all but 
SN 1998aw from the spring 1998 search were found 
at maximum light. SN 2000fr was found well 
before maximum. Thus, the new high-redshift 
dataset is likely to suffer more Malmquist bias 
than the P99 dataset. Further complicating the 
interpretation for the high-redshift supernovae is 
the fact that our new HST supernovae are spread 
over a wide range in redshift, such that a single 
brightness correction for Malmquist bias causes a 
more complicated change in the fitted cosmologi- 
cal parameters. This is unlike the situation in P99 
in which most supernovae were at z ^ 0.5. Fol- 
lowing the calculation in P99 for a high-redshift 
flux-limited SN sample we estimate that the max- 
imum Malmquist bias for the ensemble of HST su- 
pernovae is ~ 0.03 mag. However, we caution that 
it is supernovae near the flux limit which are most 
strongly biased, and therefore, that a subsample 
comprised of the highest-redshift members drawn 
from a larger flux-limited sample will be more bi- 
ased. When combined with the P99 high-redshift 
supernovae, the bias is likely to be ^ 0.02 mag 
since both samples have roughly the same statis- 
tical weight. 

As for the low-redshift SNe la, in P99 we estab- 
lished that since most of the SNe la from the H96 
flux-limited search were found near maximum, 
that sample suffered about 0.04 mag of Malmquist 
bias. On the other hand, some of the R99 SNe la 
were discovered using a galaxy-targeted technique, 
which therefore is not limited by the SN flux and 
may be more akin to a volume-limited sample (Li, 
Filippenko, & Riess 2001). Thus, the addition of 



the R99 SNe la could slightly reduce the overall 
Malmquist bias of the low-redshift sample. If we 
were to assume no Malmquist bias for the R99 
SNe la, and allowing for the fact that they con- 
tribute only ~ 1/4 the statistical weight of the 
H96 supernovae, we estimate that the Malmquist 
bias in the current low-redshift sample is roughly 
0.03 mag. 

Since Malmquist bias results in the selection of 
overly-bright supernovae at the limits of a flux- 
limited survey, and since the flux-limit can be 
strongly correlated with redshift^^, this bias can 
result in an apparent distortion of the shape of 
the Hubble diagram. This may affect estimates of 
the dark energy equation of state. The selection 
effects for the current high-redshift supernovae are 
not sufficiently well-defined to warrant a more de- 
tailed modeling of this effect than is presented 
here. However, for future work, much better con- 
trol of the selection criteria for SNe la at both 
low- and high-redshift will be required in order to 
properly estimate the impact of this small bias. 

For the current study, however, we simply note 
that since the differences in the Malmquist bi- 
ases of the high- and low-redshift subsets of SN 
are likely to be smaller in this work than in P99, 
the current results are less likely to be affected 
by Malmquist bias. Given the above estimates 
of 0.03 mag of bias in the low-redshift sample, 
and 0.02 mag of bias in the high-redshift sample, 
the difference in the biases is only 0.01 mag. To 
perform a quantitative estimate of the effects of 
Malmquist bias, we have performed a fit by ap- 
plying the mean offsets described above to each 
member of a sample in our primary subset. This 
fit is plotted in Figure 13b. The H96 supernovae 
have their magnitudes increased (made dimmer) 
by 0.04, the P99 supernovae by 0.01, and six of 
the eight HST supernovae in our primary subset 
have their magnitudes increased by 0.04. The two 
HST supernovae (SNe 1998bi, and 2000fr) which 
were found before maximum light are assumed 
not to be biased, and the other nine are offset 
by 0.04, yielding the above estimated 0.03 mag- 
nitudes for the sample. A flt with these changed 
values to the supernova peak magnitudes yields a 
flat-universe value which is different from our pri- 

'They arc 100% correlated for a single field, but this correla- 
tion can be diluted by combining fields of different depths. 
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mary fit by Qm = 0.01, and a value of Qm + 
which is different by 0.18. The best-fit value of w, 
when combined with the other cosmological mea- 
surements, is 0.03 larger. We adopt these values 
all much less than our statistical uncertainties — as 
our Malmquist bias systematic error. 

5.4. if-corrections and Supernova Colors 

The generation of the spectral template used 
for calculating /C-corrections is described in § 2.3. 
The degree to which uncertainties in the K- 
correction introduce systematic uncertainties into 
the cosmological parameters depends on whether 
or not extinction corrections are being individu- 
ally applied to supernovae. In particular, our K- 
corrections are most uncertain in the rest-frame U - 
band range of the supernova spectrum, due to lim- 
ited published spectrophotometry. As discussed 
in § 2.2, our primary fits use a spectral template 
which has a color ?7-i3=— 0.4 at the epoch of B- 
maximum. We have investigated the effects on our 
cosmology of replacing the spectral template used 
both for JC-corrections and for determining color 
excesses with a template that has U-B=— 0.5 at 
the epoch of maximum B light. 

Figure 13c shows the effect on the fitted cos- 
mology caused by using the different template for 
calculating /T-corrections when individual host- 
galaxy extinction corrections are not applied. 
These effects are very mild, indicating that our 
ii'-corrections are robust with respect to the in- 
trinsic U-B color of a supernova. Based on the 
comparison of these fits, we adopt a if-correction 
systematic uncertainty of 0.13 on Om + and of 
0.01 in w: the systematic uncertainty on the flat- 
universe value of f^M due to this effect is negligible. 

Although the effects of a different intrinsic U -B 
color on the if-corrections are mild, the effects on 
calculated color excesses are much greater. Fig- 
ure 13d shows the difference between Fit 6, where 
host-galaxy extinction corrections have been ap- 
plied using our standard color-excess values, and a 
fit where color-excess values have been determined 
assuming the intrinsic U-B color of a supernova is 
—0.5 at maximum light. As with other systemat- 
ics, the primary effect is to move the confidence 
intervals along their major axis. In this case, the 
large shift in fijvi + ^A is mainly due to the fact 
that with this bluer reference U-B color, we would 
believe that all of our z > 0.7 supernovae are suf- 



fering from an amount of host-galaxy extinction 
which is greater than that suffered by supernovae 
at lower redshift. Given that the more distant su- 
pernovae are dimmer and thus closer to our detec- 
tion limits than the moderate redshift supernovae, 
this scenario is implausible. If anything, one would 
expect the higher redshift supernovae to be less 
subject to host-galaxy extinction due to selection 
effects. Nonetheless, a value of U-B=— 0.5 at the 
epoch of _B-band maximum is currently possible 
given the [/-band information available. Only for 
those fits where extinction corrections are applied, 
we have an additional intrinsic U-B systematic er- 
ror of 0.07 on the flat- universe value of ^m, and 
a systematic error of 1.78 on Om + ^a- The sys- 
tematic uncertainty on w is 0.10. It is likely that 
these values represent an overestimate of this sys- 
tematic. 

5.5. Dust Properties 

As discussed in § 3, Phillips et al. (1999) found 
that some of the reddest supernovae at low redshift 
appear to be overcorrected for extinction given the 
standard reddening law. As shown in the lower 
panel of Figure 7, our most reddened high-redshift 
supernova (SN 1998as, which is omitted from the 
primary subset) is similarly overcorrected. One 
possible explanation is that a lower value of Rb is 
appropriate for SN la host galaxies. If we use a 
value of Rb = 3.5 (Phillips et al. 1999) rather than 
the standard value of Rb = 4.1 to perform ex- 
tinction corrections, it slightly changes the best-fit 
cosmological values for fits where extinction cor- 
rection are applied (Fit 6); this change is shown in 
Figure 13e. The best-fit value of f^M + ^^A changes 
by 0.18, and the best-fit value of w when com- 
bined with the other cosmological measurements 
changes by 0.01; this systematic has a negligible 
effect on the flat-universe value of flu- 

A related source of systematic error is possible 
evolution in the properties of the host-galaxy dust. 
To examine the scale of the effect, we consider a 
situation where dust in 2: < 0.3 spiral galaxies has 
a CardeUi, Clayton, & Mathas (1989) Ry = 3.1 
law whereas higher-redshift galaxy dust has a ra- 
tio of selective-to-total extinction that is half as 
large, i.e. Ry = 1.6. We use the Monte Carlo 
code described in Kim et al. (2003) to study the 
bias induced when an Ry = 3.1 extinction correc- 
tion is inappropriately applied to all supernovae. 
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We incorporate the rcdshift and E{B-V) distribu- 
tions of the supernovae considered in this paper 
and an E{B-V) < 0.1 cut is appHed. For an input 
cosmology of Om = 0.21 and fl^ — 0.79, we find 
a modest shift in the cosmological parameters to 
SIm = 0.25 and Oa = 0.77 without assuming a flat 
universe. 

This bias moves almost exactly along the line 
Q.M + = Ij increasing uncertainty along the 
thin axis of the error c;ontour. However, the ex- 
treme difference in dust properties considered in 
the Monte Carlo contributes a shift in the cos- 
mological parameters that is less than 1 cr of our 
quoted statistical error bars. We adopt 0.04 as the 
"dust evolution" systematic uncertainty on 
in a flat universe for those fits where host-galaxy 
extinction corrections are applied; this particular 
systematic is insignificant along the major axis of 
the confidence ellipses. 

The flat-universe value of w, when combined 
with the 2dFGRS and CMB results, increases by 
0.06 under this simple model of dust evolution. 
We adopt this as the dust evolution systematic 
on w for those fits where host-galaxy extinction 
corrections are applied. 

5.6. Gravitational Lensing 

Gravitational lensing decreases the modal 
brightness and causes increased dispersion and 
positive skewness in the Hubble diagram for high- 
redshift supernovae. The size of the effect de- 
pends on the fraction of compact objects of the 
total mass density of the universe, Om- This 
has been discussed in some detail in the litera- 
ture (Wambsganss et al. 1997; Frieman 1997; Holz 
1998; Kantowski 1998; Seljak & Holz 1999; Met- 
calf & Silk 1999; Metcalf 1999; Holz 2001; Wang, 
Holz, & Munshi 2002; Minty, Heavens, & Hawkins 
2002; Amanullah, Mortsell & Goobar 2003; Dalai 
et al. 2003; Oguri, Suto, & Turner 2003), espe- 
cially in relation to the P99 and Riess et al. (1998) 
SN datasets. A very conservative assumption of 
an "empty beam" model in a universe filled with 
compact objects allowed P99 to demonstrate that 
gravitational lensing does not alter the case for 
dark energy. 

Gravitational lensing may result in a biased de- 
termination of the cosmological parameter deter- 
mination, as discussed in Amanullah, Mortsell & 



Goobar (2003). The potential bias increases with 
the redshift of the supernovae in the sample. For 
example, for the most distant known Type la SN, 
SN1997ff at z=1.7, there is evidence for significant 
magnification, Am ~ —0.3 (Lewis & Ibata 2001; 
Mortsell, Gunnarsson & Goobar 2001; Benitez et 
al. 2002). 

As the SN sample considered in this paper does 
not reach as far, the (de)magnification distortions 
are expected to be small, in general below 0.05 
magnitudes, and less than 1% for the cases con- 
sidered in P99. To estimate the systematic uncer- 
tainties in the cosmological parameters we have 
used the SNOG package (Goobar et al. 2001) to 
simulate 100 realizations of our data sets assum- 
ing a 20% universal fraction of ^Im in compact 
objects, i.e. of the same order as the halo frac- 
tion deduced for the Milky Way from microlens- 
ing along the line of sight to the Large Magellanic 
Gloud (Alcock et al. 2000). The light beams are 
otherwise assumed to travel through space ran- 
domly filled with galaxy halos with mass density 
equally divided into SIS and NFW profiles, as de- 
scribed in Bergstrom et al. (2000). According to 
our simulations we find that (for a fiat universe) 
the fitted value of Om is systematically shifted by 
0.01 on the average, with a statistical dispersion 
ctaOm = 0.01. We adopt 0.01 as our gravitational 
lensing systematic error in the flat-universe value 
of Q.M ■ The effect on JIm + a is very small com- 
pared to other systematics, biasing the sum by 
only 0.04. 

The simulated offsets due to gravitational lens- 
ing, when combined with CMB and galaxy redshift 
distortion measurements, increase the value of w 
by 0.05; we adopt this as a gravitational lensing 
systematic on w. 

5.7. Supernova Population Drift 

In P99 we discussed in detail whether the high- 
redshift SNe la could have systematically different 
properties than low-redshift SNe la, and in partic- 
ular, whether intrinsic differences might remain af- 
ter correction for stretch. One might imagine this 
to occur if the range of the physical parameters 
controlling SN la brightnesses have little overlap 
between low- and high-redshift such that correc- 
tions applied to low-redshift are inappropriate or 
incomplete for high-redshift SNe la. Since P99, 
considerable additional work has been done to ad- 
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dress this issue. 

In addition to comparisons of stretch range 
(P99), as well as spectral (Perlmutter et al. 1998; 
Coil et al. 2000) and lightcurve (Goldhaber et al. 
2001) features, several tests performed directly 
with the P99 high-redshift SNe la have shown ex- 
cellent consistency with low-redshift SNe la. Most 
recently, in Sullivan et al. (2003) we have pre- 
sented results on the Hubble diagram of distant 
Type la supernovae from P99 that have been 
morphologically-typed with HST. We found no 
difference in the cosmological results from their 
morphologically-segregated subsamples. In par- 
ticular, E/SO galaxies — for which one expects the 
tightest possible correlation between progenitor 
mass and redshift — not only agree with the cosmo- 
logical fits using only spiral galaxies, but by them- 
selves confirm the results of P99. This is strong ev- 
idence that, while age or metallicity could in prin- 
ciple affect the brightnesses of SNe la, stretch cor- 
rection eliminates these differences. Likewise, the 
lightcurve rise-time — a possible indicator of the 
energetics of the SN explosion (see Nugent et al. 
1995; Hoeflich, Wheeler, & Thielemann 1998)— 
while initially suggested to be different between 
high- and low-redshift SNe la (Riess et al. 1999b), 
has been demonstrated to agree very well (within 
1.8 ± 1.2 days, Aldering, Knop, & Nugent 2000). 

On the theoretical side, the SN formation mod- 
els of Kobayashi et al. (1998) and Nomoto, Naka- 
mura, & Kobayashi (1999) suggest that the pro- 
genitor binary system must have [Fe/H]> —1 in 
order to produce a SN la. This would impose 
a lower limit to the metallicities of all SNe la, 
and thus limit the extent of any metallicity- 
induced brightness differences between high- and 
low-redshift SNe la. On the empirical side, the 
lack of a gradient in the intrinsic luminosities 
of SNe la with galactocentric distance, coupled 
with the fact that metallicity gradients are com- 
mon in spiral galaxies (Henry & Worthey 1999), 
lead Ivanov, Hamuy, & Pinto (2000) to suggest 
that metallicity is not a key parameter in control- 
ling SNe la brightnesses at optical wavelengths — 
though note that Lentz et al. (2000) show how it 
can affect the ultraviolet. In addition, Hamuy et 
al. (2000, 2001) find that lightcurve width is not 
dependent on host-galaxy metallicity. 

Alternatively, population age effects, including 
pre-explosion cooling undergone by the progenitor 



white dwarf and other effects linked to the mass 
of the primary exploding white dwarf have been 
suggested (for a review, see Ruiz-Lapuente 2003). 
As the local sample of SNe la represents popu- 
lations of all ages and metallicities, both effects 
can be studied locally. Several low-redshift stud- 
ies have presented data siiggcsting that SNe la in- 
trinsic luminosities (i.e., those prior to stretch cor- 
rection) may correlate with host-galaxy environ- 
ment (Hamuy et al. 1996b; Branch, Romanishin, 
& Baron 1996; Wang, Hoeflich, & Wheeler 1997; 
Hamuy et al. 2000; Ivanov, Hamuy, & Pinto 2000; 
Howell 2001; Wang et al. 2003; R99). These find- 
ings are actually encouraging, since unlike stretch 
itself, there is some hope that host-galaxy envi- 
ronment variations can be translated into physical 
parameters such as age and metallicity. These pa- 
rameters can help relate any drifts in the SNe la 
population to evolution of the host galaxies. 

More importantly for cosmology, R99 used their 
sample of 22 local SNe la to demonstrate that 
any brightness variations between SNe la in dif- 
ferent host-galaxy environments disappear after 
correction for lightcurve width. We have quan- 
tified this agreement using a larger local sample 
of supernovae compiled in Wang et al. (2003), 14 
of which have E/SO hosts and 27 of which have 
spiral hosts. We find that after lightcurve-width 
correction there can be less than a 0.01 ±0.05 mag 
offset between SNe la in local spirals and ellipti- 
cals. This indicates that lightcurve width is able 
to correct for age or other differences. 

Finally, Wang et al. (2003) demonstrate a new 
method, CMAGIC, which is able to standard- 
ize the vast majority of local SNe la to within 
0.08 mag (in contrast to ~ 0.11 mag which 
lightcurve-width corrections can attain (Phillips et 
al. 1999)). This imposes even more severe limits 
on the fraction of SNe la generated by any alter- 
nate progenitor scenario, or requires that varia- 
tions in the progenitor properties have little effect 
on whether the resulting SN can be standardized. 

The data from the new SNe la presented here do 
offer one new test for consistency between low- and 
high-redshift SNe la. The quality of our HST data 
provides measurements of the SN peak magnitudes 
and lightcurve widths rivaling those for nearby 
SNe la. This allows a direct comparison between 
the stretch-luminosity relations at low- and high- 
redshifts. Figure 14 shows that the HST high- 
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Fig. 14. — Stretch-luminosity relationship 
for low-redshift SNe (open circles) and high- 
redshift HST SNe (filled circles). Each point is 
the JC-corrected and extinction-corrected ms 
for that supernova, minus T>c, the "Hubble- 
constant-free luminosity distance" (see § 2.4), 
plotted against the stretch of that SN. The 
line drawn represents the best-fit values of a 
and M from Fit 6, the fit to all Subset 1 su- 
pernovae with host-galaxy extinction correc- 
tions applied. Note in particular that our HST 
SNe la all have low-redshift counterparts. 



Goobar, Bcrgstrom, & Mortsell 2002; Blakcslee et 
al. 2003). A direct test for extinction over a wide 
wavelength range, rest-frame B-I, have been per- 
formed by Riess et al. (2000) on a single supernova 
at 2; = 0.46, SN 1999Q, which showed no grey dust 
signature; however, see Nobili et al. (2003). Al- 
though the situation remains inconclusive, there is 
no direct evidence that "grey" dust is a dominant 
source of uncertainties. It remains an important 
issue to be addressed by future data sets including 
near-infrared observations. 

More recently, the possibility of axion-photon 
oscillations making high-redshift supernovae ap- 
pear dimmer was suggested by Csaki, Kaloper, 
& Terning (2002). This attenuation would be 
wavelength dependent, and thus could be explored 
with spectroscopic studies of high-shift sources 
(Mortsell, Bcrgstrom, & Goobar 2002). Prelim- 
inary studies of QSO spectra between z = 0.15 
and 2; = 5.3 set a very conservative upper limit on 
the possible dimming of z~0.8 supernovae to 0.2 
magnitudes (Mortsell & Goobar 2003) 

For the current data sample, the above men- 
tioned sources of systematic uncertainties are dif- 
ficult to quantify at present, but are believed to 
be subdominant in the total error budget. 



redshift supernovae are found at similar stretches 
and luminosities as the low-redshift supernovae. 
The low- and high-redshift samples are consistent 
with the same stretch-luminosity relationship, al- 
though it is primarily the low-redshift supernovae 
that require a non-zero slope for this relationship. 

5.8. Possible Additional Sources of Sys- 
tematic Uncertainties 

Other potential sources of systematic tmccr- 
tainties have been suggested. Aguirre (1999a, b) 
and Aguirre & Zoltan (2000) argued that the pres- 
ence of "grey" dust, i.e. a homogeneous inter- 
galactic component with weak differential extinc- 
tion properties over the rest-frame optical wave- 
length regime could not be rulc;(l out by the P99 
data. Since then, measurements of a SN la at 
2: ~ 1.7 (Riess et al. 2001) were claimed to rule 
out the "grey" dust scenario as a non-cosmological 
alternative explanation to the dimming of high- 
redshift supernovae; however, there remain some 
outstanding issues with this interpretation (e.g.. 



5.9. Total Identified Systematic Uncer- 
tainty 

The identified systematic errors are summa- 
rized in Table 9. Adding together these errors in 
quadrature, we obtain a total systematic error of 
0.04 on the flat-universe value of Om (along ap- 
proximately the minor axis of the confidence el- 
lipses shown in vs. plots); this is smaller 
than but approaching our statistical uncertainty 
of 0.06. The total systematic uncertainty on 
Om + f^A is 0.96 (along approximately the ma- 
jor axis of the confidence ellipses). Finally, for 
the low-extinction subset, we have a systematic 
uncertainty on constant w of 0.09, less than our 
high-side statistical uncertainty of 0.15. 

For fits with host-galaxy extinction corrections 
applied, we have to consider the additional sys- 
tematic effects of an uncertainty in the intrinsic 
value of U-B on determined color excesses, and of 
dust properties. In this case, we have a total sys- 
tematic error of 0.09 on the flat-universe value of 
or ^A, and a total systematic error of 2.0 on 
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Table 9: Identified Systematic Uncertainties 



Source of 


Systematic Uncertainty On: 


Notes 


Uncertainty 


Flat-Universe 










r^Mor JIa" 


flu + 


constant 




Fit method 


0.03 (0.5cr) 


0.80 


0.02 




Type contamination 


0.03 (0.5ct) 


0.48 


0.07 




Malmquist Bias 


0.01 (0.2(7) 


0.18 


0.03 




Intrinsic U-B: ii'-corrections 


O.OO (O.Ocr) 


0.13 


0.01 


c 


Gravitational Lensing 


0.01 (0.2cr) 


0.04 


0.05 




Systematic with host-galaxy extinction corrections: 






Intrinsic U-B: color excess 


0.07 (0.7(t) 


1.78 


0.10 


d 


Extinction Slope 


0.00 (O.Ocr) 


0.18 


0.01 


d 


Dust Evolution 


0.03 (0.3cr) 


0.02 


0.06 


d 



a: Each systematic is given as an offset from the flat- universe value of Qm, and in terms of the 
smaller side of the statistical error bar (0.06 for Fit 3 to the low-extinction subset, 0.10 for Fit 6 
to the full primary subset). 

b: This is the ofl'set on the maximum-likelihood value of w when the the fit is combined with the 
2dFGRS and CMB measurements. 

c: Only used where host-galaxy extinction corrections are not applied. 
d: Only used where host-galaxy extinction corrections are applied. 



r^M + ^^a; as discussed in § 5.4, this is likely to be 
an overestimate of the true systematic error. The 
total systematic uncertainty on constant w for the 
extinction-corrected full primary sample is 0.15. 

6. Summary and Conclusions 

1. Wc present a new, independent set of eleven 
high-redshift supernovae {z = 0.36-0.86). 
These supernovae have very high-quality 
photometry measured with WFPC2 on the 
HST. The higher quality lightcurve mea- 
surements have small enough errors on each 
E{B-V) measurement to allow an unbi- 
ased correction of host-galaxy reddening. 
We have performed improved color and K- 
corrcctions, necessary to combine WFPC2 
photometric filters with ground-based pho- 
tometric filters. 

2. The cosmological fits to flu and Oa are 
consistent with the SCP's previous results 
(P99), providing strong evidence for a cos- 
mological constant. This is a significant con- 
firmation of the results of P99 and Riess et 
al. (1998), and represents a completely new 
set of high-redshift supernovae yielding the 
same results as the earlier supernova work. 



Moreover, these results are consistent with a 
number of other cosmological measurements, 
and together with other current cosmological 
observations is pointing towards a consensus 
CIm ~ 0.3, Oa ~ 0.7 Universe. 

3. Most identified systematic errors on Om and 

$1a affect the cosmological results primarily 
by moving them along the direction where 
the statistical uncertainty is largest, that is, 
along the major axis of the confidence el- 
lipses. Systematics are much smaller along 
the minor (approximately Om — ^^a) axis 
of the confidence regions, and may be de- 
scribed by giving the systematic error on Qm 
or fl\ alone in the flat-universe case. Our 
total identified systematic error for the low- 
extinction sample analysis is 0.04 on the flat- 
universe value of Om or Oa- For flts with 
host-galaxy extinction corrections, a conser- 
vative estimate of the total identifled sys- 
tematic error is 0.09. 

In the more uncertain major axis, our total 
identified systematic error is 0.96 on ^^m+^^a 
for the low-extinction primary subset, and 
2.0 on the extinction-corrected full primary 
subset. Given the large size of these sys- 
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tematics in this direction, any conclusions 
drawn from the positions of supernova con- 
fidence ellipses along this direction should be 
approached with caution. 

4. Under the assumption of a flat universe 
with vacuum energy (constant w = — 1), 
we find a value of Om = 0.25to;o6 (sta- 
tistical) ±0.04 (identified systematic), or 
equivalently, a cosmological constant of 

= 0.75tJ5;o7 (statistical) ±0.04 (identi- 
fied systematic). This result is robust 
to host-galaxy extinction, and a fit with 
full, unbiased, individual extinction correc- 
tions applied yields a flat-universe c;osniolog- 
ical constant of flA = 0.72lQ;j^° (statistical) 
±0.09 (identified systematic). Our best con- 
fidence regions for flu versus Qa are shown 
in Figure 8. 

5. When combined with the 2dFGRS galaxy 
redshift distortion measurement and re- 
cent CMB data, we find a value for the 
dark energy equation of state parameter 
w = — l-05lo;2o marginalizing over f^M (or 
a mass density CIm = 0.27lo'o5 marginaliz- 
ing over w), under the assumptions that the 
Universe is spatially flat and that w is con- 
stant in time. The identified systematic un- 
certainty on w is 0.09. The current confi- 
dence regions on the flat-universe values of 
Om and w are shown in Figure 12. The 
supernovae data are consistent with a low- 
mass Universe dominated by vacuum energy 
{w = —1), but they are also consistent with a 
wide range of constant or time-varying dark 
energy models. 

In summary, high-redshift supernovae continue 
to be the best single tool for directly measur- 
ing the density of dark energy. This new set of 
supernovae observed with the HST confirm and 
strengthen previous supernova evidence for an ac- 
celerating universe, and show that those results 
are robust even when host-galaxy extinction is 
fully accounted for. High-redshift supernovae, to- 
gether with other cosmological measurements, are 
providing a consistent picture of a low-mass, fiat 
universe filled with dark energy. The next task for 
cosmologists is to better measure the properties of 
the dark energy, so as to further our understand- 



ing of its nature. Combinations of current cosmo- 
logical techniques have begun to provide measure- 
ments of its most general property (specifically, 
the equation of state parameter when it is assumed 
to be constant). Future work will refine these mea- 
surements, and in particular reduce the systematic 
uncertainties that will soon limit the current series 
of supernova studies. As new instruments become 
available, it will begin to be possible to relax the 
condition of a constant equation of state parame- 
ter, and to question whether the properties of the 
dark energy have been changing throughout the 
history of the Universe. 
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A. Lightcurve Data 

Tabulated below are liglitc;iirvc data for the eleven HST supernovae presented in this paper. For each event, 
there are two lightcurves, one for i?-band and one for /-band. All photometry has been color-corrected to the 
standard Bessel filters as described in § 3, using color corrections which assume the lightcurve parameters 
in Table 3. These lightcurves, together with a 7" x 7" thumbnail of the F675W WFPC2 image closest to 
maximum light, are shown in Figures 1 and 2. Note that there are correlated errors between the data points. 
For the ground-based data, there is a covariance because for a given supernova the same final reference images 
were subtracted from all other ground-based points. Similarly, the HST data include a covariance due to a 
single background model having been used for all points for a given supernova (see § 2.1). In addition to 
this, the relative photometric zeropoint magnitudes were determined separately for the ground-based and 
HST photometry; in the former case, standard stars from Landolt (1992) were used to measure magnitudes 
of secondary standard stars in the supernova field of view. In the latter case, zeropoints from Dolphin (2000) 
were used. These covariance matrices will be available from the SCP website. 

Because uncertainties are flux uncertainties rather than magnitude uncertainties, each lightcurve is pre- 
sented in arbitrary flux units. For each lightcurve, the zeropoint necessary to convert these to magnitudes 
is given. The magnitude may be calculated using the standard formula: 

m = —2.5 log/ + rrizp (Al) 

where m^p is the quoted zeropoint and / is the flux value from the table. (Because we include early-time 
and late-time lightcurve points when the supernova flux is undetected given our photometry errors, some of 
the measured fluxes scatter to negative values. Note that it is impossible to formally calculate a magnitude 
for these points, and also that flux values are the proper way to quote the data as they better reflect the 
units in which our photometry errors are approximately Gaussian.) 

The telescope used for each data point is indicated. BTC = the Big Throughput Camera on the CTIO 4m 
telescope. CTIO = the prime focus imager on the CTIO 4m telescope. WIYN the Nasmyth 2kx2k imager 
on the WIYN 3.5m telescope at Kitt Peak observatory. INT = the WFC (wide-field camera) on the INT 
2.5m telescope at La Palma. KECK = the LRIS imager on the Keck 10m telescope. NTT = the SUSI-2 
imager on the NTT 3.6m telescope at ESO. CFHT = the CFHT12K muhi-chip imager on the 3.6m CFHT 
telescope on Mauna Kea in Hawaii. Finally, HSTPC indicates data obtained from the Planetary Camera 
CCD on WFPC2. 



http: //supernova, lbl.gov/ 
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Table 11: SN1997ek-I 



Table 10: SN 1997ek-R 



Julian Dav 
-2. lUD.DUd 


Flu 


x" 


T^lo^^opc 


50780.63 


0.24 ± 1 


.27 


BTC 


50780.69 


0.57 ± 0, 


.93 


BTC 


50781.61 


-0.28 ± 1 


.05 


BTC 


50781.66 


1.22 ± 


.89 


BTC 


50781.67 


0.29 ± 0, 


.89 


BTC 


50781.72 


0.16 ± 1, 


.01 


BTC 


50810.58 


2.71 ± 1, 


.28 


BTC 


50810.59 


4.63 ± 1, 


.29 


BTC 


50810.60 


5.25 ± 1, 


.24 


BTC 


50810.67 


4.85 ± 1 


.32 


BTC 


50810.68 


5.04 ± 1 


.24 


BTC 


50810.69 


5.70 ± 1 


.28 


BTC 


50811.66 


4.34 ± 1 


.10 


BTC 


50811.68 


4.53 ± 1, 


.07 


BTC 


50811.69 


3.55 ± 1 


.22 


BTC 


50817.67 


4.92 ± 0, 


.91 


BTC 


50817.68 


5.09 ± 


.84 


BTC 


50817.69 


3.17±0 


.83 


BTC 


50817.70 


2.65 ± 


.84 


BTC 


50817.71 


3.71 ± 


.85 


BTC 


50817.72 


3.34 ± 1 


.02 


BTC 


50817.73 


4.45 ± 1 


.06 


BTC 


50817.73 


4.77 ± 1 


.04 


BTC 


50817.74 


3.10 ± 1 


.04 


BTC 


50818.92 


4.18 ± 0, 


.23 


HSTPC 


50824.77 


3.61 ± 0, 


.21 


HSTPC 


50835.67 


2.49 ± 0, 


.87 


BTC 


50835.68 


3.20 ± 


.90 


BTC 


50835.69 


2.56 ± 


.99 


BTC 


50835.70 


3.01 ± 1 


.05 


BTC 


50835.70 


3.26 ± 1 


.12 


BTC 


51165.71 


-0.05 ± 


.60 


BTC 


51165.71 


-0.67 ±0, 


.61 


BTC 


51165.74 


-0.55 ± 0, 


.71 


BTC 


51166.63 


0.44 ± 2 


.12 


BTC 


51166.65 


1.20 ± 1 


.28 


BTC 


51166.66 


-0.67 ± 1 


.49 


BTC 


51193.59 


0.47 ± 


.77 


BTC 


51193.60 


-0.86 ± 


.79 


BTC 


.")1193.61 


0.76 ± 


.70 


BTC 


51193.62 


0.18 ± 


.73 


BTC 


51194.65 


0.46 ± 0, 


.64 


BTC 



a: Zeropoint: 25.678 



Julian Day Flux^ Telescope 



-2,400,000 



OUfSiO.uU 


K ftO -1-1 AK 
O.D^ It i.40 


D L\j 


OVOl ( .oD 


Q OO -1- 1 Qn 
O.ZZ It i.oU 




OUoi (.Of 


A OT -4- 1 Qt; 
It i.OO 


£5 1 


OUoi 1 .oo 


A TCi -1-1 An 
4. / U ± i ,4U 


iS i L/ 


OUoi 1 .Oo 


^ A^ -1-1 A'i 
D.4i ± i.4o 


i) i 1^ 


oUoi ( .oy 


o.oZ ± l.ob 


BTC 


oUoi ( .bU 


4.4 f it i .DO 


BTC 


OUoi / -Di 


5. 16 ih 1.52 


BTC 


OUoi 1 .Oo 


Q ftfi -1- 1 
O.Oo It i.O^ 


i3 1 


OUoi / .04 






OUoi / .o4 


Q "31 -L 1 KCi 

o.oi ih i.oy 


BTC 


OUoi / -DO 


K (3Q -1- 1 OQ 

o.oy It i.zo 


i5 1 


oUoi ( .DO 


A QO -L 1 A A 

4.00 it i.44 


BTC 


OUoio.:.^o 


O.OO It U.iD 


ilO i 


t; nc 1 O 7/1 


o m -1- 1 Tn 


W i 1 IN 


c nc 1 o Td 
OUoiy. 1 D 


o.Uo It l.Do 


VV i 1 IN 


oUoiy. io 


A 1 O -L 1 on 

4. io it i.yu 


W i I IN 


OUoiy. li) 


1 71 -1- 1 ftn 

i. ( i It l.DU 


W 1 I IN 


ouoiy.oi 


4:.Oi It 1.0 ( 


W 1 I IN 


50819.82 


3.84 lb 2.09 


WIYN 


50824.78 


3.89 ±0.16 


HSTPC 


50835.72 


2.72 lb 1.96 


BTC 


50835.73 


3.06 ± 2.05 


BTC 


50846.74 


1.54 ± 0.09 


HSTPC 


50858.84 


0.75 ± 0.07 


HSTPC 


50871.95 


0.46 ± 0.06 


HSTPC 


51072.07 


0.50 lb 0.57 


KECK 


51072.07 


0.35 lb 0.58 


KECK 


51072.07 


0.69 lb 0.58 


KECK 


51072.11 


0.31 lb 0.55 


KECK 


51072.11 


0.94 lb 0.58 


KECK 


51072.12 


-0.23 ± 0.57 


KECK 


51101.99 


-0.37 ± 0.54 


KECK 


51102.00 


0.51 ± 0.58 


KECK 


51102.00 


0.58 ± 0.59 


KECK 


51102.05 


1.20 lb 0.75 


KECK 


51102.06 


1.53 lb 0.90 


KECK 


51126.93 


-0.04 lb 0.06 


HSTPC 


51134.26 


0.06 ± 0.05 


HSTPC 


51165.70 


-0.66 ± 1.15 


BTC 


51165.72 


0.21 ± 1.06 


BTC 


51165.73 


-0.44 ± 1.12 


BTC 


51193.64 


0.01 ± 1.12 


BTC 


51193.65 


-0.28 ± 1.13 


BTC 


51193.67 


-0.46 lb 1.50 


BTC 


51194.59 


0.99 lb 1.17 


BTC 


51194.60 


1.34 lb 1.30 


BTC 


51194.60 


0.73 lb 1.15 


BTC 



a: Zeropoint: 24.801 
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Table 12: SN 1997eq-R 



Julian Day Flux" Telescope 



-2.400,000 



50780 


.60 


0, 


.01 


± 


0, 


12 


BTC 


50780 


.66 


0, 


.21 


± 


0, 


12 


BTC 


50781 


.60 


-0, 


.08 


± 


0. 


,10 


BTC 


50781, 


.63 


0, 


.19 


± 


0. 


.10 


BTC 


50781, 


.68 


0, 


.09 


± 


0. 


.10 


BTC 


50781, 


.72 


0, 


.14 


± 


0. 


.11 


BTC 


50810, 


.61 


1, 


.76 


± 


0. 


,12 


BTC 


50810 


.62 


1 


.80 


± 


0, 


,12 


BTC 


50810 


.63 


1 


.88 


± 


0, 


13 


BTC 


50810 


.64 


1 


.87 


± 


0, 


,11 


BTC 


50810 


.70 


1 


.91 


± 


0, 


,12 


BTC 


50810, 


.71 


1 


.82 


± 


0, 


,11 


BTC 


50811, 


.70 


1, 


.78 


± 


0. 


.10 


BTC 


50818 


.34 


2, 


.23 


± 


0, 


,28 


INT 


50818 


.36 


1 


,98 


± 


0, 


24 


INT 


.■50810 


.85 


1 


,69 


± 


0, 


,05 


HSTPC 


50821 


.66 


2 


.14 


± 


0, 


,54 


WIYN 


50821 


.67 


1 


.79 


± 


0, 


,39 


WIYN 


50835 


.41 


0, 


.85 


± 


0, 


,13 


INT 


50835 


.42 


0, 


.87 


± 


0, 


,18 


INT 


50835 


.43 


0, 


.85 


± 


0, 


,34 


INT 


50843, 


.68 


0, 


.37 


± 


0. 


.18 


WIYN 


50843, 


.70 


0, 


.02 


± 


0. 


,40 


WIYN 


50846 


.81 


0, 


.32 


± 


0, 


,02 


HSTPC 


50855 


.82 


0, 


.18 


± 


0, 


,02 


HSTPC 


50863 


.82 


0, 


.12 


± 


0, 


,02 


HSTPC 


51165 


.56 


0, 


,01 


± 


0, 


12 


BTC 


51165 


.61 


0, 


,01 


± 


0, 


,41 


BTC 


.■51165 


.62 


-0, 


,61 


± 


0, 


67 


BTC 


51165 


.64 


0, 


.00 


± 


0, 


,12 


BTC 


51193, 


.58 


-0, 


.03 


± 


0. 


.10 


BTC 


51193, 


.63 


0, 


.02 


± 


0. 


.09 


BTC 



a: Zeropoint: 23.284 



Table 13: SN 1997eq-I 



Julian Day 
-2,400,000 




Flux" 


Telescope 


50818.37 


1, 


.15 ±0. 


.50 


INT 


50818.38 


1, 


.05 ±0. 


32 


INT 


50818.39 


1, 


.20 ± 0. 


,32 


INT 


50818.41 


0, 


.94 ± 0. 


,49 


INT 


50818.43 


1 


.20 ± 0, 


,48 


INT 


50818.46 


1 


.05 ± 0, 


,25 


INT 


50819.87 


0, 


.91 ± 0, 


,03 


HSTPC 


50821.68 


0, 


.93 ± 0. 


.35 


WIYN 


50821.69 


0, 


.83 ±0. 


.41 


WIYN 


50821.70 


0, 


.65 ±0. 


.38 


WIYN 


50824.90 


0, 


.86 ±0, 


,02 


HSTPC 


50835.54 


0, 


.59 ± 0, 


,27 


INT 


50835.56 


0, 


,13 ± 0, 


,29 


INT 


50835.58 


-0, 


.11 ± 0, 


,50 


INT 


50846.82 


0, 


.38 ± 0, 


,02 


HSTPC 


50855.83 


0, 


.27 ±0, 


,02 


HSTPC 


50863.83 


0, 


.22 ± 0. 


.01 


HSTPC 


51165.57 


0, 


.03 ± 0. 


,29 


BTC 


51165.60 


0, 


.06 ± 0, 


,34 


BTC 


51165.63 


0, 


.07 ± 0, 


,20 


BTC 


51165.65 


0, 


.06 ± 0. 


.17 


BTC 


51193.58 


-0, 


.10 ± 0. 


.17 


BTC 



a: Zeropoint: 22.388 



Tabic 14: SN 1997ez-R 



Julian Day 
-2,400,000 


Flux" 


Telescope 


50780.75 


—0.41 ± 1 


.15 


BTC 


50780.82 


—0.88 ± 


.96 


BTC 


50781.74 


-1.46 ± 1, 


.01 


BTC 


50781.79 


0.29 ± 1, 


.18 


BTC 


50781.79 


1.09 ±0 


.96 


BTC 


50811.77 


6.05 ± 1, 


.04 


BTC 


50811.77 


3.90 ± 1 


.89 


WIYN 


50811.77 


5.82 ± 1 


.03 


BTC 


50811.78 


5.62 ± 1 


.02 


BTC 


50811.78 


5.82 ± 2 


.22 


WIYN 


50811.79 


3.97 ±4 


.73 


WIYN 


50811.81 


5.97 ± 1, 


.04 


BTC 


50811.81 


4.83 ± 1, 


.16 


BTC 


50817.84 


5.51 ± 1, 


.22 


BTC 


50817.85 


7.72 ± 1 


.63 


BTC 


50817.86 


4.58 ± 2 


.15 


BTC 


50818.70 


4.93 ± 1 


.13 


INT 


50818.72 


5.04 ± 1 


.09 


INT 


50819.06 


4.96 ± 


.25 


HSTPC 


50824.97 


3.65 ± 


.22 


HSTPC 


50835.66 


4.69 ± 1, 


.49 


INT 


50835.67 


2.88 ± 1, 


.68 


INT 


50835.81 


1.81 ± 1, 


.49 


BTC 


50835.82 


-0.07 ± 1, 


.66 


BTC 


50835.83 


0.52 ± 1 


.70 


BTC 


51193.75 


-0.14 ± 


.74 


BTC 


51193.76 


0.37 ± 


.69 


BTC 


51193.76 


0.00 ± 1 


.08 


BTC 


51193.77 


-1.23 ± 


.85 


BTC 


51193.78 


-0.20 ± 0, 


.83 


BTC 


51193.79 


-0.21 ± 0, 


.78 


BTC 


51193.80 


-1.80 ± 1 


.63 


WIYN 


51195.73 


-1.37 ± 1 


.26 


WIYN 


51195.75 


-0.21 ± 1 


.40 


WIYN 


51195.77 


-0.58 ± 1 


.18 


WIYN 


51195.78 


-0.92 ± 1, 


.36 


WIYN 



a: Zeropoint: 25.688 
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Tabic 15: SN 1997ez-I 



Julian Day 
-2,400,000 


Flux" 


Telescope 


50816.74 


2.05 ± 1 


.90 


BTC 


50816.76 


4.83 ± 2, 


.03 


BTC 


50816.77 


4.64 ± 1 


.89 


BTC 


50816.78 


6.11 ± 1, 


.90 


BTC 


50816.78 


5.02 ± 2, 


.02 


BTC 


50816.85 


6.84 ±2, 


.14 


BTC 


50818.63 


4.19 ± 2 


.23 


INT 


50818.65 


4.24 ± 1 


.55 


INT 


50818.66 


4.12 ± 1 


.54 


INT 


50818.68 


4.30 ± 1 


.54 


INT 


50819.07 


5.23 ± 0, 


.18 


HSTPC 


50820.79 


4.42 ± 1 


.56 


WIYN 


50820.81 


5.69 ± 1, 


.50 


WIYN 


50820.83 


3.92 ± 1 


.46 


WIYN 


50820.84 


4.22 ± 1, 


.42 


WIYN 


50820.86 


6.08 ± 1 


.67 


WIYN 


50820.87 


3.26 ± 1 


.70 


WIYN 




4.07 ± 


.17 


xlo 1 r^O 


50835.60 


5.27 ± 1 


.77 


INT 


50835.61 


0.53 ± 2 


.03 


INT 


50835.63 


5.55 ± 1, 


.94 


INT 


50835.64 


5.62 ± 2, 


.52 


INT 


50835.84 


3.39 ± 2 


.13 


BTC 


50835.85 


1.78 ± 2, 


.23 


BTC 


50835.86 


-0.47 ± 2 


.56 


BTC 


50846.55 


1.77 ± 


.09 


HSTPC 


50858.98 


1.00 ± 


.08 


HSTPC 


50871.89 


0.48 ± 


.04 


HSTPC 


51189.97 


0.80 ± 1 


.13 


WIYN 


51189.98 


-0.74 ± 1, 


.22 


WIYN 


51190.00 


-0.20 ± 1, 


.35 


WIYN 


51191.90 


-0.54 ± 1 


.34 


WIYN 


51191.92 


-1.64 ± 1 


.16 


WIYN 


51191.93 


0.15 ± 1 


.28 


WIYN 


51194.70 


-3.19 ± 2 


.44 


BTC 


51194.71 


-1.06 ± 2 


.73 


BTC 


51194.72 


-0.60 ± 2 


.43 


BTC 


51194.73 


-0.52 ± 2 


.81 


BTC 


51194.74 


-1.26 ± 2, 


.28 


BTC 


51194.75 


-0.84 ± 2, 


.49 


BTC 


51194.76 


-0.27 ± 1 


.90 


BTC 


51194.77 


-2.00 ± 2 


.19 


BTC 


51194.78 


-1.89 ± 2 


.02 


BTC 


51194.78 


-1.58 ± 2 


.61 


BTC 


51194.79 


-0.68 ± 2 


.38 


BTC 



a: Zeropoint: 24.954 



Table 16: SN 1998as-R 



Julinn Da\' Flux" 'I'rloropc 

-2.1II(J.()II(I 



50872.63 


— 


, 10 ± 0, 


10 


BTC 


008 I A.OO 


— 


,07 i 0, 


09 


BTC 


008 I A.O I 





.06 ± 0, 


09 


BTC 


50872.72 


— 


.07 ± 0, 


10 


BTC 


50872.73 


-0, 


.06 ± 0, 


,11 


BTC 


50873.57 


0, 


.06 ±0, 


,11 


BTC 


50873.58 


0, 


.03 ±0, 


,10 


BTC 


50895.58 


2, 


.33 ±0, 


,12 


BTC 


50895.62 


2, 


.47 ±0, 


,15 


BTC 


50896.58 


2, 


.64 ± 0, 


,12 


BTC 


50899.70 


2 


.24 ± 0, 


,12 


BTC 


50904.68 


2, 


.15 ± 0, 


,11 


BTC 


50904.69 


2, 


.05 ± 0, 


,10 


BTC 


50904.70 


2, 


.20 ± 0. 


,10 


BTC 


50904.71 


1, 


.95 ± 0, 


,11 


BTC 


50904.72 


2, 


.00 ± 0, 


,10 


BTC 


50912.29 


1, 


.42 ± 0, 


,04 


HSTPC 


50935.01 


0, 


.33 ± 0, 


,02 


HSTPC 


50948.52 


0, 


.25 ± 0, 


,02 


HSTPC 


50963.17 


0, 


.19 ± 0, 


,02 


HSTPC 


51193.83 


0, 


.06 ± 0, 


,08 


BTC 


51193.84 


-0, 


,07 ± 0, 


,08 


BTC 


51193.86 


0, 


,04 ± 0, 


,08 


BTC 


51196.03 


0, 


.21 ± 0, 


,13 


WIYN 


51196.04 


-0, 


,19 ± 0, 


,12 


WIYN 


51196.05 


-0 


,11 ± 0, 


,16 


WIYN 



a: Zeropoint: 23.139 



Tablo 17: SM998as-I 

Julian Day Flux'^ Telescope 



-2,400,000 



50912.31 


9, 


,24 ± 


,21 


HSTPC 


50924.07 


7, 


,27 ± 


,19 


HSTPC 


50932.65 


1 


,95 ± 1 


,56 


WIYN 


50935.02 


4, 


,86 ± 


,17 


HSTPC 


50948.53 


2 


,57 ± 


,14 


HSTPC 


50963.19 


1 


,79 ± 


,12 


HSTPC 


51194.86 


-1 


,02 ± 


.98 


BTC 


51194.87 


0, 


,60 ± 1, 


.12 


BTC 


51196.93 


-0, 


,55 ± 1, 


.23 


WIYN 


51196.94 


0, 


,73 ± 1 


.12 


WIYN 


51196.96 


-1 


,44 ± 1 


.28 


WIYN 


51280.50 


0, 


,53 ± 1 


.60 


BTC 


51280.51 


-2, 


,08 ± 1 


,50 


BTC 


51280.51 


0, 


,67 ± 1 


.50 


BTC 


51280.52 


0, 


,60 ± 1, 


.33 


BTC 


51280.53 


1, 


.32 ± 1 


.45 


BTC 


51280.54 


0, 


,72 ± 1, 


.46 


BTC 



a: Zeropoint: 24.788 
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Tabic 18: SN 1998aw-R 



Julian Day 
-2,400,000 




Flux" 


Telescope 


50513.71 





.08 ± 0, 


14 


BTC 


50513.73 


— 


.08 ± 0, 


,16 


BTC 


50513.75 


0, 


.06 ±0, 


,13 


BTC 


50514.71 


0. 


.08 ± 0, 


,14 


BTC 


50517.74 


-0, 


.19 ±0, 


,14 


BTC 


50517.76 


0. 


.04 ± 0, 


,16 


BTC 


50518.79 





.31 ± 0, 


,17 


BTC 


50518.81 


— 


.02 ± 0, 


,17 


BTC 


50872.56 


-0 


.03 ± 0, 


,21 


BTC 


50872.59 


-0 


.03 ± 0, 


,22 


BTC 


50873.73 


— 


.03 ± 0. 


,18 


BTC 


50873.74 


— 0, 


.09 ±0, 


,15 


BTC 


50895.60 


0, 


.02 ± 0, 


,16 


BTC 


50895.64 


0, 


.55 ± 0, 


,16 


BTC 


50896.58 


0. 


.67 ±0, 


,15 


BTC 


50896.60 





.39 ± 0, 


,16 


BTC 


50899.69 





.89 ± 0, 


,15 


BTC 


ouyu4.Do 


1 


.87 ± 0, 


,14 




50904.64 


1 


.66 ± 0, 


,14 


BTC 


50904.65 


1 


.75 ± 0, 


,13 


BTC 


50904.66 


1, 


.82 ± 0, 


,14 


BTC 


50904.67 


1, 


.82 ± 0, 


,14 


BTC 


50912.03 


2, 


.53 ± 0, 


,07 


HSTPC 


50922.11 


2, 


.11 ± 0, 


,06 


HSTPC 


50927.56 


2, 


.05 ± 0, 


,38 


BTC 


50927.57 


1, 


.80 ± 0, 


,34 


BTC 


50927.60 


1 


.69 ± 0, 


,36 


BTC 


50927.61 


0, 


.96 ± 0, 


,41 


BTC 


50929.64 


1 


.48 ± 0. 


,28 


WIYN 


50929.65 


1, 


.06 ± 0, 


,33 


WIYN 


50929.67 


1 


.90 ± 0, 


,31 


WIYN 


509.33.07 


1 


,32 ± 0, 


04 


HSTPC 


.■50947.71 





,58 ± 0, 


,03 


HSTPC 


50961.83 


0, 


,30 ± 0, 


,03 


HSTPC 


51192.96 


-0, 


,19 ± 0, 


,26 


WIYN 


51192.98 


-0, 


,14 ± 0, 


39 


WIYN 


51193.00 


0, 


,18 ± 0, 


28 


WIYN 


51193.02 


-0, 


,14 ± 0, 


,24 


WIYN 


51193.03 


-0, 


,29 ± 0, 


,28 


WIYN 


51279.60 


0, 


,01 ± 0, 


,13 


BTC 


51279.61 


0, 


,04 ± 0, 


,14 


BTC 


51279.63 


-0, 


,04 ± 0, 


,12 


BTC 


51279.66 


0, 


,01 ± 0, 


,13 


BTC 


51280.56 


0, 


,14 ± 0, 


,16 


BTC 


51280.57 


0, 


,17 ±0, 


,15 


BTC 



a: Zeropoint: 23.536 



Table 19: SN 1998aw-I 



Julian Day 
-2,400,000 




Flux" 


Telescope 


50513.76 


-0, 


,33 ± 0, 


,25 


BTC 


50514.74 


-0, 


,10 ±0, 


,22 


BTC 


50514.76 


-0, 


,12 ± 0, 


,21 


BTC 


50514.78 


0, 


,06 ± 0, 


23 


BTC 


50518.73 


0, 


,18 ± 0, 


,42 


BTC 


50518.75 


-0, 


,08 ± 0, 


,34 


BTC 


50912.04 


1 


,79 ± 0, 


,05 


HSTPC 


50922.12 


1, 


,67 ±0, 


,05 


HSTPC 


50929.70 


1, 


,50 ±0, 


,49 


WIYN 


50930.71 


1, 


,80 ±0, 


,46 


WIYN 


50933.08 


1 


,23 ±0, 


,03 


HSTPC 


50947.73 


0, 


,80 ± 0, 


,03 


HSTPC 


50961.84 


0, 


,53 ± 0, 


,03 


HSTPC 


51194.03 


-0, 


,07 ± 0, 


,32 


WIYN 


51194.05 


-0, 


,26 ± 0, 


,51 


WIYN 


51195.97 


-0, 


,21 ± 0, 


,32 


WIYN 


51195.98 


0, 


,13 ± 0, 


,27 


WIYN 


51196.00 


0, 


,10 ± 0, 


,29 


WIYN 


51196.02 


0, 


,05 ± 0, 


,27 


WIYN 


51279.59 


-0, 


,03 ± 0, 


,21 


BTC 


51279.62 


-0, 


,06 ± 0, 


,25 


BTC 


51279.64 


0, 


,15 ± 0, 


,21 


BTC 


51279.65 


0, 


,01 ± 0, 


,23 


BTC 


51279.66 


0, 


,19 ± 0, 


,25 


BTC 


51280.55 


0, 


,14 ± 0, 


,31 


BTC 


51280.57 


-0, 


,02 ± 0, 


,28 


BTC 


51280.59 


-0 


,30 ± 0, 


,29 


BTC 


51280.60 


0, 


,09 ± 0, 


,29 


BTC 



a: Zeropoint: 22.874 
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Tal)l(> 20: SX l<)<)8ax-R 



Julian Day 




Telescope 


-2,400,000 







50138.65 


-0, 


.03 ± 0, 


,09 


CTIO 


50138.67 


-0, 


.09 ± 0, 


10 


CTIO 


50159.64 


-0, 


.09 ± 0. 


,08 


CTIO 


50159.66 


0, 


.03 ± 0, 


,07 


CTIO 


50160.67 


0, 


.01 ±0, 


,07 


CTIO 


50160.68 


0, 


.02 ±0, 


,06 


CTIO 


50168.59 


-0, 


.03 ±0, 


,07 


CTIO 


50168.65 


0, 


.14 ± 0, 


,06 


CTIO 


50169.64 


0, 


.13 ± 0, 


,15 


CTIO 


50169.67 


-0, 


,01 ± 0, 


,08 


CTIO 


50432.83 


-0, 


,06 ± 0, 


06 


CTIO 


50453.84 


-0, 


.01 ± 0, 


,08 


CTIO 


50454.77 


0, 


.01 ±0, 


,06 


CTIO 


50459.82 


-0, 


.02 ± 0, 


,04 


CTIO 


50459.83 


-0, 


.02 ± 0, 


,05 


CTIO 


50459.84 


0, 


.02 ± 0, 


,05 


CTIO 


50490.79 


0, 


,01 ± 0, 


,06 


BTC 


50490.79 


0, 


,07 ± 0, 


06 


BTC 


50490.80 


-0, 


,04 ± 0, 


,06 


BTC 


50490.80 


-0, 


,04 ± 0, 


,06 


BTC 


50513.71 


-0, 


,03 ± 0, 


,06 


BTC 


50514.72 


-0, 


.06 ± 0, 


,06 


BTC 


50872.54 


0, 


,72 ± 0, 


,12 


BTC 


50872.57 


0, 


,58 ± 0, 


,12 


BTC 


50873.53 


0, 


,84 ± 0, 


,17 


BTC 


50873.55 


0, 


,95 ± 0, 


,10 


BTC 


50895.52 


1 


,42 ± 0, 


,09 


BTC 


50895.55 


1 


,06 ± 0, 


,19 


BTC 


50895.71 


1 


,24 ± 0, 


,07 


BTC 


50896.53 


1 


,14 ± 0, 


,10 


BTC 


50900.70 


1, 


,14 ±0, 


,07 


BTC 


50900.71 


1, 


.04 ± 0, 


,07 


BTC 


50904.59 


0, 


,91 ± 0, 


,06 


BTC 


50904.60 


0, 


,84 ± 0, 


,06 


BTC 


50904.61 


0, 


,81 ± 0, 


,06 


BTC 


50904.62 


0, 


,84 ± 0, 


,06 


BTC 


50904.63 


0, 


,89 ± 0, 


,06 


BTC 


50911.96 


0, 


,59 ± 0, 


,03 


HSTPC 


50922.04 


0, 


,31 ± 0, 


,02 


HSTPC 


50933.00 


0, 


,18 ± 0, 


,02 


HSTPC 


50947.65 


0, 


,09 ± 0, 


,01 


HSTPC 


50961.23 


0, 


,09 ± 0, 


,01 


HSTPC 


51193.80 


-0, 


,00 ± 0, 


,05 


BTC 


51193.81 


-0 


,00 ± 0, 


,05 


BTC 


51193.82 


-0, 


,01 ± 0, 


06 


BTC 


51279.52 


-0 


,01 ± 0, 


,08 


BTC 


51279.57 


0, 


,11 ± 0, 


,08 


BTC 


51280.61 


0, 


,06 ± 0, 


,06 


BTC 



a: Zeropoint: 22.922 



Table 21: SN1998ax-I 



Julian Day 

-2,400,000 



Flux" Telescope 



50911.97 


1.95 ± 0.10 


HSTPC 


50922.05 


1.62 ± 0.10 


HSTPC 


50933.01 


1.18 ± 0.06 


HSTPC 


50947.66 


0.75 ± 0.05 


HSTPC 


50961.24 


0.47 ± 0.04 


HSTPC 



Tabic 22: SN1998ay-R 



Julian Day 


Flux" 


Telescope 


-2,400,000 







50521.85 





,02 ± 0, 


50 


WIYN 


50521.86 





,17 ± 0, 


56 


WIYN 


50872.54 


2 


,11 ± 1, 


,08 


BTC 


50872.57 


1, 


,27 ±0, 


,97 


BTC 


50873.53 


0, 


,57 ± 1, 


,81 


BTC 


50873.55 


-0 


,70 ± 1, 


04 


BTC 


50895.52 


5 


,69 ± 0, 


90 


BTC 


50895.55 


6, 


,69 ± 1, 


91 


BTC 


50895.71 


6, 


,10 ± 0, 


,78 


BTC 


50896.53 


6, 


,70 ± 1, 


,24 


BTC 


50900.70 


5 


,74 ± 0, 


,76 


BTC 


50900.71 


6, 


,74 ± 0, 


,91 


BTC 


50904.59 


5, 


,48 ± 0, 


,78 


BTC 


50904.60 


5, 


,64 ± 0, 


,75 


BTC 


50904.61 


5 


,61 ± 0, 


,78 


BTC 


50904.62 


5 


,76 ± 0, 


,82 


BTC 


50904.63 


5 


,91 ± 0, 


,79 


BTC 


50912.16 


3, 


,11 ± 0, 


,20 


HSTPC 


50923.99 


1 


,58 ± 0, 


17 


HSTPC 


51193.80 


-0, 


::- 


60 


BTC 


51193.81 


0, 


,61 ± 0, 


,48 


BTC 


51193.82 


0, 


,53 ± 0, 


,64 


BTC 



o: Zeropoint: 25.093 



Table 23: SN1998ay-I 



Julian Day 
-2,400,000 



Flux" Telescope 



50912.17 


1.56 ±0.08 


HSTPC 


50924.00 


0.96 ±0.07 


HSTPC 


50934.68 


0.61 ± 0.04 


HSTPC 


50948.59 


0.40 ± 0.04 


HSTPC 


50967.81 


0.26 ± 0.04 


HSTPC 


o: Zeropoint: 


23.685 





Zeropoint: 23.685 
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Tabic 26: SN 1998be-R 











Julian Day 




Flux"* 


Telescope 










-2,400,000 


















50490, 


,86 


0, 


,49 ± 0, 


,55 


BTC 


Table 24: SN 1998ba-R 




50490, 


,87 


— 0, 


,39 ±0, 


,54 


BTC 


tT 1 1 1 K 1 11 (I \^ 






T^lp^^OI5C 


50513, 


,83 


-0, 


,02 ± 0, 


52 


BTC 








50513, 


,84 





,15 ± 0, 


,.54 


BTC 


oyjo 1 o. 1 iJ 


Q 


Q'l -l. n no 


BTC 


50514, 


,83 





,53 ± 0, 


,60 


BTC 


ouo 1 o.oy> 


Q 


no + n no 


BTC 


,■50514, 


,86 


— 


,51 ± 0, 


,53 


BTC 


iJVjO 1 O.O X 


Q 


ni + n no 


BTC 


50517, 


,88 





,33 ± 0, 


,70 


BTC 




Q 


fiQ + n no 


BTC 


50517, 


,90 


— 


,26 ±0, 


,71 


BTC 


iJUO I o.oo 


Q 


ni -1- n ns 


BTC 


50517, 


,90 


0, 


,69 ±0, 


,81 


BTC 


OKJO I O-Ci 


— 


03-1-0 09 


BTC 


50518, 


,86 


0, 


,22 ± 0, 


,62 


BTC 


ouoyo. 1 




50 i 14 


BTC 


50518, 


,87 


0, 


,57 ±0, 


,66 


BTC 






.64 lb 0.15 


BTC 


50872, 


,74 


— 


,75 ± 0, 


,91 


BTC 


50899.75 


1, 


.52 ±0.11 


BTC 


50872, 


,89 


1 


,36 ± 0, 


,93 


BTC 


50899.84 


1 


.43 ± 0.14 


BTC 


50873, 


,87 





,63 ± 0, 


,53 


BTC 


50899.90 


1 


.20 ± 0.21 


BTC 


50895, 


78 


4 


,22 ± 0, 


,70 


BTC 


50900.74 


1 


.54 ± 0.10 


BTC 


50895, 


,84 


5 


,34 ± 0, 


,88 


BTC 


50900.75 


1 


.32 ± 0.10 


BTC 


50899, 


,75 


7, 


,13 ± 0, 


,79 


BTC 


50904.77 


1 


.36 ± 0.11 


BTC 


50899, 


,82 


6, 


,98 ± 0, 


,91 


BTC 


50904.78 


1, 


.20 ± 0.11 


BTC 


50900, 


,76 


4, 


,64 ± 0, 


,65 


BTC 


50904.79 


1, 


.42 ± 0.13 


BTC 


50904, 


,73 


6 


,58 ± 0, 


,65 


BTC 


50904.80 


1, 


.30 ± 0.09 


BTC 


50904, 


,74 


6, 


,90 ± 0, 


,67 


BTC 


50904.81 


1, 


.34 ±0.11 


BTC 


50904, 


,75 


6, 


,31 ± 0, 


,72 


BTC 


50912.10 


0, 


.79 ± 0.03 


HSTPC 


50904, 


,75 


7, 


,32 ± 0, 


,73 


BTC 


50923.12 


0, 


.41 ± 0.02 


HSTPC 


50904, 


,76 


8, 


,29 ± 0, 


,76 


BTC 


50933.21 


0, 


.22 ± 0.02 


HSTPC 


50904, 


,86 


7, 


,95 ± 0, 


,89 


BTC 


50947.12 


0, 


.12 ± 0.01 


HSTPC 


50912, 


,23 


5, 


,73 ± 0, 


,25 


HSTPC 


50961.90 


0, 


.12 ± 0.01 


HSTPC 


50923, 


,19 


2, 


,11 ± 0, 


,18 


Xlk3 1 X^kv 


51258.01 


-0, 


.15 ± 0.11 


WIYN 


50932, 


,74 


2, 


,04 ± 0, 


,89 


WIYN 


51279.82 


0, 


.07 ±0.08 


BTC 


50932, 


,77 


1, 


,38 ± 0, 


,93 


WIYN 


51279.85 


-0, 


.05 ± 0.10 


BTC 


50934, 


,08 


0, 


,73 ± 0, 


,12 


HSTPC 


51280.69 


-0, 


.02 ± 0.07 


BTC 


50949, 


,00 


0, 


,76 ± 0, 


,13 


HSTPC 


51280.70 


0, 


.03 ± 0.06 


BTC 


50962, 


,17 


0, 


,21 ± 0, 


,13 


HSTPC 


a: Zeropoint: 22.779 




51279, 


,68 


-0, 


,16 ± 0, 


,67 


BTC 










51279, 


,71 


0, 


,31 ± 0, 


,68 


BTC 










51279, 


,75 


0, 


,21 ± 0, 


,73 


BTC 










51279, 


,77 


-0, 


,30 ± 0, 


,79 


BTC 










a: Zeropoint: 25.350 














Table 27: SN 1998bc-I 












Julian Day 




Flux" 


Telescope 










-2,400,000 


















50514, 


,85 


-0, 


,21 ± 0, 


,83 


BTC 


Tm]i1(> 9.-,- in()X]iM-T 




50514, 


,87 


-1 


,02 ± 0, 


,78 


BTC 










50518, 


,84 


2, 


,00 ± 0, 


,90 


BTC 


Julian Day 




Flux" 


Telescope 


50518, 


,85 


1, 


,47 ± 0, 


,86 


BTC 


-2,400,000 








50518, 


,85 


0, 


,31 ±0, 


,82 


BTC 


50907.82 


3, 


.18 ± 1.99 


W i 1 IN 


50912, 


,25 


3, 


,66 ± 0, 


,18 


HSTPC 


50907.83 


3, 


.96 ± 1.75 


VV i 1 IN 


50923, 


,20 


2 


,19 ± 0, 


,17 


HSTPC 


50907.84 


6, 


.80 ± 1.81 


W i 1 IN 


50932, 


,80 


2 


,35 ± 1, 


,10 


WIYN 


50907.85 


6, 


.04 ± 2.36 


WIYN 


50932, 


,85 


2 


,26 ± 0, 


,92 


WIYN 


50912.11 


5, 


.74 ±0.22 


HSTPC 


50934, 


,09 


1 


,13 ± 0, 


,09 


HSTPC 


50923.13 


3, 


.95 ±0.21 


HSTPC 


50949, 


,01 


0, 


,80 ±0, 


,08 


HSTPC 


50933.22 


2, 


.81 ±0.12 


HSTPC 


50962, 


,19 


0, 


,37 ±0, 


,08 


HSTPC 


50947.13 


1, 


.57 ±0.10 


HSTPC 


51279, 


,69 


0, 


,81 ± 0, 


,89 


BTC 


50961.92 


1, 


.37 ± 0.10 


HSTPC 


51279, 


,70 


0, 


,49 ± 0, 


,87 


BTC 


51279.83 


-1 


.51 ± 1.00 


BTC 


51279, 


,72 


1 


,51 ± 0, 


,73 


BTC 


51279.84 


0, 


.88 ± 1.09 


BTC 


51279, 


,73 


-0, 


,02 ± 0, 


,71 


BTC 


51280.69 


-1 


.04 ± 0.83 


BTC 


51279, 


,76 


0, 


,62 ± 0, 


,83 


BTC 


51280.71 


0, 


.66 ±0.72 


BTC 


51279, 


,77 


0, 


,58 ± 0, 


,85 


BTC 


51280.72 


-0, 


.06 ± 0.68 


BTC 


51280, 


,64 


-0, 


,87 ±0, 


,82 


BTC 


51280.73 


0, 


.13 ± 0.68 


BTC 


51280, 


,64 


0, 


,36 ± 0, 


,84 


BTC 


a: Zeropoint: 24.477 




51280, 


,65 


0, 


,12 ± 0, 


,73 


BTC 










51280, 


,66 


-0, 


,13 ± 0, 


,78 


BTC 










51280, 


,67 


1, 


,24 ± 0, 


,76 


BTC 










51280, 


,68 


-0 


,62 ± 0, 


,76 


BTC 



a: Zeropoint: 24.384 
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Table 30: SN2000fr-R 



Tabic 28: SN 1998bi-R 

Julian Day Flux" Telescope 



-2,400,000 



50138.79 


-1.04 ±0, 


,91 


CTIO 


50138.82 


0.85 ±0, 


,86 


CTIO 


50168.80 


-0.68 ±0, 


66 


CTIO 


50490.86 


0.40 ± 0, 


,49 


BTC 


50490.87 


-0.09 ± 0, 


,48 


BTC 


50513.83 


0.26 ± 0, 


,51 


BTC 


50513.84 


-0.10 ± 0, 


,53 


BTC 


50514.83 


-1.06 ± 0. 


,58 


BTC 


50514.86 


-0.05 ±0, 


.50 


BTC 


50517.88 


0.13 ±0, 


.65 


BTC 


50517.89 


-0.11 ±0, 


.60 


BTC 


50517.89 


0.93 ±0, 


,60 


BTC 


50517.90 


-0.29 ± 0, 


,68 


BTC 


50517.90 


-0.35 ± 0, 


,74 


BTC 


50872.89 


0.52 ± 0, 


,81 


BTC 


50873.87 


0.60 ± 0, 


,51 


BTC 


50895.78 


3.15 ± 0. 


,63 


BTC 


50895.84 


3.11 ± 0. 


.79 


BTC 


50899.75 


4.93 ± 0. 


.65 


BTC 


50899.82 


4.28 ± 0, 


.70 


BTC 


50900.76 


4.44 ± 0, 


,55 


BTC 


50904.73 


6.10 ± 0, 


,61 


BTC 


50904.75 


5.30 ± 0, 


,61 


BTC 


50904.75 


5.38 ± 0, 


,64 


BTC 


50904.76 


6.21 ± 0, 


,66 


BTC 


50904.86 


5.27 ±0, 


,77 


BTC 


50910.15 


5.27 ±0, 


.20 


HSTPC 


50922.18 


3.75 ± 0, 


.18 


HSTPC 


51279.71 


0.94 ±0. 


.73 


BTC 


51279.74 


0.63 ± 0, 


.67 


BTC 


51279.75 


-1.14 ± 0, 


.68 


BTC 


51279.77 


0.47 ±0, 


.76 


BTC 


a: Zeropoint: 


25.213 







Julian Day Flux'^ Telescope 



-2,400,000 



51671.77 


1 


,02 ± 0, 


,07 


KECK 


51671.77 


1 


.05 ± 0, 


,07 


KECK 


51671.78 


1, 


.06 ±0, 


,07 


KECK 


51671.78 


0, 


.99 ±0, 


,07 


KECK 


51679.98 


1, 


.66 ±0, 


,04 


HSTPC 


51692.91 


1 


,43 ± 0, 


,03 


HSTPC 


51706.26 


0, 


,73 ± 0, 


,02 


HSTPC 


51718.04 


0, 


,39 ± 0, 


,01 


HSTPC 


51733.86 


0, 


,16 ± 0, 


,01 


HSTPC 


52014.72 


-0, 


,01 ± 0, 


,07 


NTT 


52014.73 


-0, 


.08 ± 0, 


.07 


NTT 


52014.74 


0, 


.04 ± 0, 


.08 


NTT 


52014.75 


-0, 


.04 ± 0, 


.06 


NTT 


52014.76 


-0, 


.04 ± 0, 


,07 


NTT 


52014.77 


-0, 


,08 ± 0, 


,10 


NTT 


52014.78 


-0, 


,07 ± 0, 


,09 


NTT 


52014.79 


-0, 


,04 ± 0, 


,10 


NTT 


52014.80 


-0, 


,16 ± 0, 


,14 


NTT 


52376.98 


0, 


,01 ± 0, 


,04 


CFHT 


52376.99 


-0. 


.00 ± 0, 


.03 


CFHT 


52377.04 


0, 


.01 ± 0, 


.04 


CFHT 


52377.05 


-0, 


.02 ± 0, 


.04 


CFHT 


52382.01 


0, 


.03 ± 0, 


.05 


CFHT 


52384.98 


-0, 


,00 ± 0, 


.09 


CFHT 


52386.85 


-0, 


,14 ± 0, 


.10 


CFHT 



a: Zeropoint: 22.998 



Tabic 29: SN 1998bi-I 



Julian Day 
-2,400,000 



Flux" Telescope 



50910.16 


2.07 ±0.06 


HSTPC 


50922.20 


1.83 ±0.06 


HSTPC 


50931.99 


1.25 ± 0.04 


HSTPC 


50946.38 


0.54 ± 0.03 


HSTPC 


50966.88 


0.20 ± 0.02 


HSTPC 


a: Zeropoint: 


23.685 





Ta]Ac 31: SX 2()()()fr-I 



Julian Day 


Flux" 


Telescope 


-2,400,000 







51641.99 


0, 


,03 ± 


,04 


CFHT 


51664.95 


0, 


,40 ± 


,05 


CFHT 


51664.99 


0, 


,40 ± 


,06 


CFHT 


51672.86 


1 


,14 ± 


,02 


HSTPC 


51679.97 


1 


,59 ± 


.03 


HSTPC 


51692.91 


1, 


.46 ± 0, 


.03 


HSTPC 


51706.20 


1, 


.02 ±0, 


.03 


HSTPC 


51717.98 


0, 


,66 ± 


.02 


HSTPC 


51733.79 


0, 


,40 ± 


.02 


HSTPC 


51997.93 


0, 


,05 ± 


,06 


CFHT 


51997.94 


0, 


,01 ± 


,06 


CFHT 


51997.99 


0, 


,19 ± 


.05 


CFHT 


51998.00 


0, 


.03 ± 0, 


.06 


CFHT 


51998.01 


0, 


.08 ± 0, 


.06 


CFHT 


52376.96 


0, 


.04 ± 0, 


.06 


CFHT 


52376.97 


-0, 


.06 ± 0, 


.06 


CFHT 


52377.00 


0, 


,13 ± 


.06 


CFHT 


52377.00 


-0, 


,09 ± 


.06 


CFHT 


52377.01 


-0, 


,01 ± 0, 


.06 


CFHT 


52377.03 


0, 


,01 ± 0, 


.07 


CFHT 



a: Zeropoint: 22.805 
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